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INTRODUCTION 


The present work comprises basically three essays, one for each 
observational mode of Geometric Satellite Geodesy. 

The first one, published separately as No. 185 in this Report series, 
"On the Geometric Analysis and Adjustment for Optical Satellite Observa- 
tions," is about the optical observations mode. In that, the geometrical 
analysis of the problem yielded a well-known regression model for the 
adjustment of the observations along with a suitable and convenient 
metric for the least-squares optimum criterion. In the present treatment 
of the problem, considerable emphasis has been placed on the geometry, 
and no additional element has been introduced for the adjustment. 

The second essay, which is in this Report (Section 1), is about the 
deterrainantal loci (critical configurations) for range networks. This 
problem has been treated before in (Rinner [1966]) and (Blaha [1971]), 
the latter being an extension of the former along similar guidelines. 

The present treatment approaches the topic from the viewpoint of 
analysis. An attempt is made to found the results theoretically by using 
elements of the theory of invariants. To the writer’s opinion the present 
treatment is not the most elegant one. Some partial results from an attempt 
for a pure geometrical approach, which remained unfinished, hint that a 
geometrical approach would be more elegant. 

The third essay, also in this Report (Section 2), is about the deter- 
minantal loci for networks of stations observing range differences (e.g., 
by means of geoceivers). In order to find out these determinantal loci, 
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an analytical way is followed, similar to that used for the case of range 
networks. Although the problem for "continuously" measured range differences 
presents some similarities in treatment with that of ranges, it is different 
from that in many aspects. 

A number of appendices is included which constitute a substantial 
part of the whole work. 

As for the mathematical language used, it can be said that formal 
mathematical language has been avoided. 
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1. DETERMINANTAL LOCI FOR RANGE NETWORKS 
1 . 1 Introductory Remarks 
1.1.1 Recovery of the Satellite-Position Points 

In the range observations mode each participating station Pi at an 
event [Ej , Qj|tj] observes the length of the distance (P^Qj), i.e. 

r i j = (PiQj) • 

This distance determines a spherical surface with center at P± and radi- 
us r±j, which may be considered as a locus of the satellite-position 
point Qj.. Thus Qj is determined as the intersection of spherical sur- 
faces with centers at the station-position points P± and radii 
the lengths of the distances of Qj from P±. In order to formulate it 
analytically, let [x±, y±, zi] be the Cartesian coordinates of Pi and 
[Xj, Yj, Zj] of Qj, with respect to a coordinate system tied to the 
solid earth. Then for the determination of Qj , one has the system of 
equations of the form: 

(Xj-x i ) 2 + (Yj-y i ) 2 + (Zj-z i ) 2 -rij = 0, (±=1,2,...) (1) 

where [x±, yi , z±] are considered to be known here, and the range of the 
index i remains to be determined for the unique determination of 
Qj (Xj , Yj , Zj). From the implicit function theorem [see Appendix A] it 
is inferred that a system of three equations of the form (1) may be 
solved for Xj , Yj, Zj in terms of x±, y±, z±, (i*l, 2, 3) provided cer- 
tain conditions are satisfied. Thus if the system. 
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0 


(Xj-xi) 2 +(Yj-y 1 ) 2 +(Zj-z 1 ) 2 -rij = 

(X j -X2) 2 +(Y j -y 2 ) 2+ (Zj-Z2) 2 -r2j = 0 (2) 

( x j -X3) 2 + (Yj -y3) 2 + (Zj -Z3) -rjj = 0 
does have at least one solution-point* and 

Xj-xi Yj -yi Zj-*! 

X j“ x 2 Y j'Y2 Z j" z 2 (3) 

X j -x 3 Y j-y3 2 j“ z 3 

then it can be solved for Xj , Yj , Zj in terms of Xi , y^, r^j, 

(i=l, 2, 3), in the neighborhood of a solution-point, (a point which 
annuls (2). 

Although the implicit function theorem assures the feasi- 
bility of the solution of the system (2) with respect to X j , Y j , Zj , it 
does not explain the set of the solution-points. This is in fact 

an algebraic geometry problem, and one may take the occasion to say that 
the whole deal about range and range-dif ference networks in the present 
work belongs to the jurisdiction of algebraic geometry. 

The set of solution-points of (2) is treated in Appendix B, where 
it is demonstrated that under condition (B15) a set of points which 
annul (2) exists, and under (3), that set contains two points symmetric 
with respect to the plane of the centers. Thus in order to tie the satellite- 
position points to the solid earth, three known stations observing simul- 
taneously are necessary and sufficient to do that. Therefore, in a series 

* It is noted that here and throughout this work the field of concern 
is that of real numbers. 
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of range observations events where three stations observe simul- 
taneously, one has no redundant* information at all. For redun- 
dant information at least four stations observing simultaneously 
are necessary, provided their configuration is not a degenerate 
one. When more than four stations observe simultaneously, 
they are sufficient but not necessary for redundant infor- 
mation from the observations. 

1.1.2 Introduction of the Coordinate System 

Suppose that there are m stations, , ( i = 1,2, . .., n ) 

observing simultaneously a series of range observations events 
[Ej, Qj | tj], (j = 1, 2 , . * , m ) . The configuration of the points 
Pj^, Qj is a polyhedron, not necessarily convex, in the extended 
field of geometric geodesy. The polyhedron must be "built" through 
the observations. One may consider that polyhedron being con- 
structed by piecing together tetrahedra, as building blocks. 

To start constructing the polyhedron, a basis must be taken upon which the 
which the rest is founded. It is desirable for the basis to be 
characterized as necessary and sufficient for this purpose. Such 
a basis is the 3-dimensional simplex , that is, the tetrahedron. 


* To put it formally, redundant information from a set of satellite 
observations for geometric satellite geodesy 3 is that which is beyond the 
necessary and sufficient for the unique determination of the satellite 
position points, (extension points of the field of geometric geodesy). 
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The coordinate system to be introduced, is attached to the tetrahedron 
that serves as basis. This is usually done by considering a coordinate 
system and fixing certain coordinates in such a way that a non-degener- 
ate tetrahedron is determined. If for an example, one constrains the 
three coordinates of P 1# two of the coordinates of P2 and one of the 
coordinates of P3, i.e.. 



the fixed tetrahedron is formed from the points (see Figure 1 ); 

0 (0,0,0) - origin of the coordinate system, 

P 1 Cxi, yi, zi) _ t he point, the coordinates of which are all fixed, 
C p 2)y z ^2* z 2^ ‘ t * ie projection on the yz-coordinate plane of 

the point P2 whose coordinates y and z are 
fixed, 

(P3) z ( 0 , 0 , Z3) - the projection on the z-axis of the point P3 

whose coordinate z is fixed. 

Some people choose xj = yj = z\ =■ yj ~ Z 2 = z 3 = 0 . For numerical cal- 
culations this choice is preferable; however, in a theoretical treatment 
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it is not advisable , for it destroys the symmetry of the expressions, 
and consequently deprives one from a great advantage for elegant mani- 
pulations. Therefore in the following treatment, the constrained coor- 
dinates will be left unspecified as long as there is no need to do so. 
1,1*3 Minimum Number of Events for a Unique Solution . 

Consider now the configuration of the points 
Pi, (i = 1, 2, ..., n), Q j > (j - 1, 2, m), and let 

Pi = ( x i Yi> z i)» i = n 

Qj - (Xj j Yj , Z j ) , j = 1, 2, m, 

be the analytic representations of the points P^, Q-j with respect to a 
cartesian coordinate system, which has been introduced through an admis- 
sible set of constraints, say those quoted above, i.e., 

*1. yi» Zl» Y2, *2, 23 = constants. 

Suppose that except the constrained coordinates, all the others are 
unknown. Thus there are (3n-6) unknown station-coordinates. One re- 
calls that n>4 and then (3n-6)>6. A question of interest in the follow- 
ing discussion is the minimum number of events for a unique determina- 
tion of the unknown station -coordinates . At each event three unknown 
quantities are introduced, namely the coordinates of the satellite- 
position point, while n equations result from the simultaneous observa- 
tions of the n stations, which are of the form (1). It was proved above 
that three equations are, (under certain conditions), necessary and suf- 
ficient for the determination of the satellite-position point. There- 
fore at each event (n-3) equations are redundant. For a unique deter- 
mination of the 3n-6 unknown station-coordinates, one has to find out 
min m under the condition 
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That is 


m(n-3)>3n-6, ns4. 


(4) 


min m = 


if 3 {Sif} is an inte « er 

[3.( n .~2)] + i; if 3-£. n .l i~ X is not an integer, 
(n-3) (n-3) 


(5) 


where stands for the integral part of 3 (n-2)/ (n-3) . 

(n-3) 

Although n may be at least theoretically, any positive integer >4, 
this is not the case for m, for the sequence q = 3 (n-2) /(n-3) decreases 
monotonicly, converging to 3 as n->-oo, i.e., 

li m 3 (n-2) _ 
n-^oo (n-3) 

Thus one has 

min m = 6 for n = 4 

= 5 " n = 5 (6) 

=4 " n > 5,(i.e. n = 6, 7, 8, ...). 

1.1.4 Patterns of Observations 

By pattern of observation, it is meant the way a set of stations 
observes a set of satellite-position points. For a given number of sta- 
tions, each observational pattern is a combination of some basic observa- 
tional patterns. In the following, a formula will be set up, giving the 
number of the basic observational patterns of any number of stations. 
Although such a formula does not seem to possess any practical interest, 
trying to obtain it one gains an insight to the situation. 

Let there exist n stations. It was shown above that the minimum 
number of stations participating in an event must be four. Therefore the 

n stations may observe either in groups of 4, or 5, , or (n-1) or all 

together. The total number of basic observational patterns then is 
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N = C 4 + C 5 + 


+ C 


n 

n -1 


+ 


C 


n 

n 


= 2 n - (C q + cj + + C§) 

= 2 n -(l + n + n ( n - ] ) 4. n (n~l) 
2 6 


N = 2 n -(n+l) •(l+ n ^- n ~ 1 ^ ) 

6 

The last formula gives 

N = 1 for n = 4 

= 6 " n = 5 

= 22 " n = 6 

= 64 " n = 7 

= 163 " n = 8 


( 7 ) 


( 8 ) 


1 . 2 Definition of Determinants 1 Loci 
Consider a set of n stations Pi(z ~ 1 , 2, .... n) observing in some 
pattern m events, corresponding to the satellite-position points 
Qj, (j = 1, 2, .... m) . Think of the satellite-position point Qj being 
observed simultaneously by nj stations, where nj>4. Then for each Qj 
there exist nj observation equations of the form ( 1 ), i.e. 
fjil = ( x j - x i!) 2+ CYj-y ll ) 2 +C z j-z il ) 2 - r i 1 j = 0, 


fj^2 ~ j 2 ^ + 2 } " r -*-23 "" 


(9) 


(x r xi ■p'Pj-n n ) 2 *< z j- z i n ) 2 -n j - 0, 

J J J J J 

where ij, i 2 » . i n ^ are indices from the set {1, 2, 3, n) which 

correspond to the stations observing simultaneously Qj, and 

j = 1, 2, 3, ..., m. Totally there are 
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( 10 ) 


> ) - ri|+n2 +n 3 + ' • + ‘ ' * +n m 

observation equations, and 

y = (3n+3m)-6 = 3(n+m-2) (11) 

unknowns, namely Xj , Yj, Zj, (j=l,...,m) and xj, y it z if (i-l,..,n) with 
six of them constrained. It is supposed that v>y; otherwise the system 
is under-determined. Consider the system of y equations out of v. From 
the implicit function theorem one has that if a solution point to this 
system exists*, then it can be solved uniquely for the y unknown vari- 
ables in a neighborhood of that solution-point, provided the jacobian 
of the functions of the left-hand parts of this system does not vanish 
in that neighborhood, (see Appendix A, relation A3) . Since one may 
choose y equations out of v in as many ways as Cy = v!/y!(v-y)!, in or- 
der to include all these cases, the above condition for the jacobian is 
stated as follows. The jacobian matrix of the functions of the left- 
hand parts of the system (9) (v equations in y unknowns) is of rank y. 


i.e. , 


/ 9 (fjii, •• .fjin-0 \ 

rank | - 1 = rank JM=y=3(n+m-2) 

\3(X j ,Y j ,Z j ,x i ,y i ,z i )/ 


( 12 ) 


i—l,2,..., n 

j=l,2, . . .,m 

If the rank of the above jacobian matrix is less than y, one cannot 
find y independent functions out of v, that means that additional rela- 
tions exist between the variables. Since the rank of the jacobian matrix 
is an invariant under linear transformations (see Appendix C), the addi- 
tional relations between the variables have to be invariants too , under 


*0ne recalls that in the implicit function theorem (see Appendix A) the 
existence of a solution-point is assumed. 
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linear transformations, and therefore they represent geometric relations 
(see [Tsimis, 1972] section 2. 2.2. 3). The totality of those additional 
relations, (for which the jacobian matrix is of rank less than p) con- 
stitute the determinantal loci .* Therefore in order to find out the 
determinantal loci, one has to examine when the jacobian matrix is of 
rank less than p. Although the jacobian matrix is huge, (vxp), each of 
its rows has at most six non-zero elements out of p, and when one takes 
the determinants of the submatrices pxp to check for the rank, pivotal 
condensation techniques reduce those determinants considerably (see 
Appendices D and E) . 


1.3 Determinantal Loci and Adjustment of the Observations 
Consider again the system (9) being represented by one of its equa- 
tions, in order to simplify the notation, i.e., 

fji = (X j -x i ) 2 + (Y r y i ) 2 + (Z ;j -z i ) 2 - ri ? =0 (13) 

where 


jell, 2, ... ,m) 

ie{ij ,i 2 » • *■ > i n ^ ^ ( 1 > Z , 3 , . . . , n } * 

By considering rjj as unknowns in addition to Xj, Yj 
take the differential form at (13). Then 


(14) 


'J ’ 


: i> Yi> 


'i » 


2 (Xj -xi) (dXj -dxj) + 2 (Yj -yi) (dYj -d yi ) + 2 (Zj -z ± ) (dZj -dz ± ) = Zr^dr^ 


(15) 


or 


(Xj-xj) 

r ij 


(dXj-dxi) 




(Yj-yj] 

r ij 


(dYj -dy-jj 


+ (Zj-zi) 

r ij 


(dZj -dz-^) -dri j ^ 


(16) 


*The name is justified from the fact that each of these additional rela- 
tions annuls the determinants of all pxp submatrices of the jacobian 
matrix. 


- 11 - 


. t - \rO .,0 -O O O O OO s i j i 

Refer this differential form to the point (Xj ,Yj ,Zj jX-j^yi ,Zi,rjLj) and let 
o b o b 

drij = rij-rij = (rij-rij)-(rij-rij) = vij-lij 
where r-n is the true value, r^j is the observed one v-y = r^-rij and 


ij 
o b 

lij = rij-rij. 

Then (16) yields 
,o o 




.o o 


O o. 


'^i^dX r dx k ).iIi^(dY r dy 1 )*15i^[d2 r d2 i ) = Vjj-ljj, (17) 

r ij r ij r ij 


11 

or in matrix form 

L = -AX+V 


(18) 


where 

v^l = (lij ) j v V l = ( v ij)> (p+6)^l = (dXj ,dYj ,dZj ,dxi ,dyi,dzi) 
and results from the jacobian matrix Jm (see relation [12]), by 

dividing each of its rows by the corresponding r^j . It is noted that 
since rij/O the rank of does not change by that operation and conse- 

w 

quently 

rank A = rank . (19) 

Under the assumptions E(V) = 0 and Cov(V) = E(VV^) = I where is 
unknown, relation (18) is a linear statistical model. From the theory of 
the linear model one has that the best (minimum-variance) linear (linear 
functions of lij) unbiased estimate of X is given by least squares and 
is unique iff rank A =* p . By recalling relations (19) and (12) one con- 
cludes that a unique least squares solution exists iff 
rank A = rank = p, 

i.e., if there exist no determinantal loci. That is why the determin- 
antal loci are studied in the present work. 
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1 . 4 Determinantal Loci of Specified Observational Patterns 

1.4.1 General Remarks 

In this section specified observational patterns will be examined 
in order to find out their determinantal loci. That is, given the num- 
ber of stations and the way they observe (observational pattern), the 
corresponding jacobian matrix will be searched for the cases when its rank 
is less than the number of the unknowns* Despite the fact that there 
exist a large number of cases to be examined, a systematic procedure may 
cover them all, at least by enumeration, for they constitute a denumer- 
able set. 

It is recalled [see section 1.1.1] that in general , four ground 
stations are necessary and sufficient for each event. This must especi- 
ally be stressed , for it is the core of the whole problem. Thus the 
case of four ground stations is treated first. No specification of the 
coordinates which are to be constrained is made until a general formula 
is reached that is susceptible to all the specifications. The theory of 
determinants is used extensively; however, explicit expressions 
of the determinants are given instead of symbolic, sacrificing compact- 
ness for an easier understanding of the manipulations involved. 

1.4.2 The Case of Four Ground Stations 

Let Pj = (x 1? y l3 z 1 D,-P 2 = fx 2 , y 2 . p 3 H ( x 3> ^3’ z 3^’ 

P4 = (X4, 74, X4) be four simultaneously observing ground stations, and 
Qj = (Xj, Yj, Zj), (j = 1, 2, . m) be a set of satellite-position 
points observed by those stations. Then the functional relationships 
on which the statistical model is based, have the form 
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( 20 ) 


fij s (X j -x i ) 2 +(Y j -y i ) 2 + (Z : j-z i ) 2 = r?j 
i » 1, 2, 3, 4 
j = 1, 2, . . . , m, m>6 
[see section 1.2 and 1.3). 

This system contains 

v = 4m>24 equations 


and 


( 21 ) 


3m+3x4 = 3(m+4) unknown variables, six of which are to be con- 
strained, so that there actually are 

p=3 (m+2) unknown variables. (22) 

Nevertheless, while taking the jacobian matrix (J M ) for the system (20), 
no unknown variable is considered constrained. Thus the jacobian matrix 
of the system (20) is of dimension 
4mx3(m+4) . 

According to what was said in section 1.3, the least squares solution is 
unique iff 


rank J M = p = 3 (m+2) (23) 

while for 

rank < y = 3 (m+2) (24) 

one has the determinantal loci of the configuration. Therefore in order 
to find out the determinantal loci, one has to consider the determinants 
of all yxy submatrices of the jacobian matrix J^, where the value of 
y=3(m+2) is determined from the value of m which is obtained by equating 
the number of unknowns to the number of equations, i.e. 



4m = 3 (m+2) => m=6. 

(25) 

Thus 

y = 3 (m+2) = 24 

(26) 
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It should be remembered that m=6 it the minimum number of events for a 


unique solution in the case of four ground stations, (see section 1.1.3, 
tion [6]). Manipulating with determinants of that order (i.e. 24-th) 
seems to be an appalling work; however they may be reduced in order by 
pivotal condensation in a systematic way, which leads to a sort of con- 
densation at the jacobian matrix, being called here condensed Jacobian 
matrix. In order to avoid interruption of the main objective by annoy- 
ing details, the derivation of the condensed jacobian matrix is given 
in Appendix E, from where the final form is recalled in Table 1. 


X 1 

yi 

z i 

x 2 

y 2 

z 2 

+ 

(Xj-x^a} 

+ 

( Y r y i )A i 

+ 

(Zj-Zi)Aj 

(X r x 2 )4 

(Yj-y 2 )Aj 

(Z j -z 2 )A^ 


x 3 

>'3 

z 3 

x 4 1 

1 

y 4 

z 4 

+ 

i 

j (X r x 3 )A^ 

+ 

(Yj-y 3 )4 

+ 

(Z r z 3 )ij 

1 

(Xj-x 4 )Aj 

(Yj-y 4 )Aj 

(Zj-z 4 ) &| 


Table 1. The Condensed Jacobian Matrix for Four Stations 


In Table 1: 


X r x 2 Yj-y 2 Z r Z 2 


x j Y j z j 1 

x 2 y 2 z 2 1 

Xj _x 3 Y j-y3 Z j" z 3 

— 

x 3 Y3 z 3 1 

Xj-x 4 Yj -y 4 Z r z 4 


*4 y 4 Z 4 1 

Xj-xi Yj -y i Zj-zj 


Xj Yj Zj 1 

x i yi z i 1 

Xj-x 3 Yj -y 3 Zj-z 3 


x 3 X3 z 3 1 

Xj-x 4 Yj -y 4 Zj-z 4 


x 4 y 4 z 4 1 


(j = l,2,...,m), (27) 
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rela- 




(j = l,2, . . . ,m) 


(27) 

(cont.) 


1 
1 
1 
1 

1 
1 
1 
1 

The above expression of the condensed jacobian matrix is a compact symme- 
tric form, where the constrained coordinates are not specified. Each time 
a set of admissible constraints is selected, one has to delete the corres- 
ponding columns from the above expression. Because of the symmetry it is 
very flexible to specifications. The determinants of the 24x24 submatri- 
ces of the original jacobian matrix, which are taken for the determinantal 
loci , are equal to the determinants of the corresponding 6x6 submatrices 
of the condensed jacobian matrix. Therefore for the determinantal loci 
one has to consider the determinants of all 6x6 submatrices of the con- 
densed jacobian matrix. To fix ideas, suppose that the coordinates 
X 1 >Yl > z l>/2 > z 2» z 3 are constrained. Deleting the corresponding columns 
from the condensed jacobian matrix, it becomes an mx6 array, i.e.. 


4 - 


x r x i Vn Z r z i 

x r x 2 v r y 2 Zj-z 2 
X j ~ x 4 Yj-y 4 Zj z 4 


X j 

X 1 


Y j 

y \ 


z i 


x 2 Y 2 z 2 

x 4 Y 4 z 4 


aJ 


Xj-xj Yj -yj Zj-zj 
Xj-x 2 Yj -y 2 Zj -z 2 
Xj-x 3 Yj-y 3 Zj-z 3 


X j 

X 1 

x 2 

x 3 


>1 

Y2 

Y3 


(compare Appendix E, relations [E14]). 


Z i 

Z 1 

z 2 

z 3 


x 2 

x 3 

y 3 

x 4 

Y4 

z 4 

-(Xj-x 2 )4 

+ (Xj-x 3 )A| 

+ ( Y j -Y3)4 

-(Xj-x 4 )Aj 

-(Yj-y 4 )Aj 

-(Zj-z 4 )Aj 


Ta ble 2. The Condensed Jacobian Matrix for Four Stations When the 
coordinates Xj, yj, z\, y 2 , z 2 , z 3 are constrained 

For the determinantal loci one must examine the zeroes of the determin- 
ants of the 6x6 submatrices. For a representative of these determinants. 
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consider that of the first six rows, i.e., 

-(X r x 2 )A2 (Xj-XjJAj (Y r y 3 )A3 -(Xj-x^aJ -(X r y 4 )Aj -(Zj-z 4 )aJ 

-(X2~X2)A2 (X2-x 3 )A 3 (Y 2 -y 3 )A 3 -(X 2 -x 4 )A 4 -(Y2-y 4 )A 4 -(Z 2 -z 4 )A 4 

-(X 3 -X2)A2 (X3-X3) A 3 (Y 3 -y 3 )A 3 -(X 3 -x 4 )A 4 -(Y 3 -y 4 )A 4 -(Z 3 -z 4 )A 4 

A | y ^ 6 = 

-(X 4 -x 2 )A 2 (X 4 -x 3 )A 3 (Y 4 -y 3 )A 3 -(X 4 -x 4 )A 4 -(Y 4 -y 4 )A 4 -(Z 4 -z 4 )A 4 

-(X 5 -x 2 )A2 (X 5 -x 3 )A 3 (Y 5 -y 3 )A 3 -(X 5 -x 4 )A 4 -(Y 5 -y 4 )A 4 -(Z 5 -z 4 )A 4 

-(X 6 -x 2 )A2 CX 6 -x 3 )&3 (Y 6 -v 3 )A 3 -(X 6 -x 4 )A 4 -(Y 6 -y 4 )A 4 -(Z 6 -z 4 )A 4 
(Check with the expression [E 17 ] of Appendix E) . 

In terms of fundamental invariants (see Appendix F) each of these 


determinants may be expressed as follows: 




,£< 

••° 6 ' D l(P 2 ) < 


l) c 1 (- 1 ) 2 


*2 11 

V 1 


iX 32 T JJ 1 

"p Y j3 1 


1*4 \ H 1 

■ X j4 'V Z i ’* * 

f35 Z J5 1 

iV V Z j6 




( 29 ) 


This expression is obtained by two imbedded steps of Laplace expansion, 
namely: 

1 st step-expansion according to the last three columns, and 

2 nd step-expansion of the cofactor of each term in the first step 
according to its last two columns. 

01 and 02 result from the two imbedded Laplace expansions respec- 
tively. The summation is a double one. The outer summation being taken 
over the pj selections of row indices taken three at a time, and the 
inner one over the p 2 selections of the three complementary row indices 
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for each of the p 1 selections taken two at a time. Thus summation (29) 
extends over 

PjP2 = C 3 C 2 = 3V3T 2T5T = 60 terms 
These determinants are relative invariants {see Appendix C) and 
when equated to zero, a hypersurface results which is an absolute in- 
variant , for multiplication by a factor would have no significance in 
its equation (see Appendix F). One observes that each of these deter- 
minants is a functional of twenty-four variables. Since all the vari- 
ables are coordinates of points in the geometric space, one may analyze, 
so to speak, that hypersurface into 3-dimensional loci. Each of these 
loci is a part of the hypersurface and all together constitute the hy- 
persurface. To carry out that analysis a systematic procedure which 
fits the geometric characteristics of the problem will be beneficial. 
Thus consider all the variables involved falling in two groups: (A) 

the coordinates of the satellite-position points, and (B) the coordin- 
ates of the station-position points. Then two categories of loci will 
be considered: (1) the loci for the satellite-position points, and (2) 

the loci for the station-position points. 

(1) Determinantal Loci for the Satellite-Position Points 
Consider the equation 


I A I j 1 "* ° ( ' 30) 

where jj, is a se l ec ti° n of six indices from (1,2,. 

If (30) is satisfied for any six satellite-position points, then the 
rank of the jacobian matrix (table 3) is less than six and consequently 
one has a determinantal locus the analytic representation of which is 
(30). This is a hypersurface in 24-dimensional space. However, one 
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may obtain 3 -dimensional loci on this hyper surf ace by keeping the coor- 
dinates of all the points but one constant. One observes from (29) 
that the left-hand side of (30) is a functional of 

i) second degree with respect to the coordinates Xj , Yj , Zj of 
any of the satellite-position points, 

ii) sixth degree with respect to the coordinate X2 f 

iii) fifth degree with respect to the coordinates x 3 , y 3 , and 

iv) fourth degree with respect to the coordinates X 4 , y 4 , 24 - 
Thus the three dimensional loci on that hypersurface for different 
points are generally different. Here the three dimensional loci on the 
hypersurface (30), for a satellite-position point will be considered. 
Consider, without loss of generality, the hypersurface 



with all the coordinates involved being kept constant, except the coor- 
dinates of one of the satellite-position points, say X 6 , Y 6 , Z 6 . Then 
l A ll 6 =0 represents a second degree surface as a locus of the point 

Q = (X 6 , Y 6 , Z 6 ). This surface passes through the four station-posi- 
tion points and the first five satellite-position points Qj, Q 2 , Q 3 , Q4, 
Q 5 . As a matter of fact these nine points lie on that surface, for by 
plugging the coordinates of any of them in the place of X^, Y^, Z^, the 
last row of |&li .,6 vanishes. Thus if for any selection 
out of (1,2,... ,m} the equation (30) is satisfied, which is equivalent 
to saying that all the satellite-position points lie on the second 
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degree surface defined by the four station-position points and any five 
satellite-position points, then the rank of the jacobian matrix is less 
than six. Thus, the four station-position points and any five satel- 
lite-position points determine through (30) a determinantal locus in 
the three dimensional space (a second degree surface) for the rest 
points . 

Now one proceeds to find out the linear (first degree) three dimen- 
sional determinantal loci. 

Consider the expansion (29) of va ^ ue 

determinant is zero if each of its expansion-terms vanishes. By equat- 
ing to zero each factor of the typical expansion-term in (29) one ob- 


tains the following determinantal loci. 


(i) By putting 


1 X9 1 

X- 1 ~ X 2' x i = 0, U = 1,2, .. . ,m) , it is con- 

1 1 1 


eluded that the plane perpendicular to the x-axis, passing through 


P 2 = (X 2 > y 2 > z 2 ) a determinantal locus for the satellite-position 

points. In other words if all of them lie in that plane, all 
| a| j j > . . . j vanish, and then the rank of the jacobian matrix is less 
than six. 


(ii) By equating the next factor to zero, i.e., 
x 3 X3 1 

X. Y-: 1=0 (31) 

J 2 J 2 

X- Y- 1 ( j , J any two indices from {l,2,...,m}) 

, J 3 23 

one has the col linearity condition of the projections on the (x,y)- 
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Thus if all the satellite-position points lie in a plane through 

P 3 e (x 3 , y 3 , z 3 ) and perpendicular to the (x,y)-plane, then relation 

(311 is satisfied for any - - from (1 ,2, . . . ,m} , and consequently each 

J2>^3 

plane of the pencil of planes with axis the projection line of P 3 to 
the (x,y)-plane, is a determinantal locus for the satellite-position 
points, (see Figure 1). 

(iii) Equating the third factor to zero one has 


x 4 

X4 

z 4 

l 


x. 

Y- 

Z. 

l 


^4 

H 

u 


= 0 (32) 

X Ss 

V 

Y- 

z. 

l 


Js 

^6 

Z h 

l 

(for any from { 1, 2, . . . ,m}j 


This is a condition for four points in 3-dimensional space to be line- 
arly dependent, that is for the points. 
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P 4 E (x 4 »y 4 > z 4 ) » Qj 4 = CXj 4 > Y : j 4 ,z j4 ), Qj ^ = (x j 5 ’ Y j 5 *Zj 5 ^ Qj 6 = 

to be coplanar . Therefore if all the satellite-position points lie in 
a plane through the point P 4 2 (x 4 ,y 4 ,z 4 ) relation (32) is satisfied 


for any triple •? ^ from l,2,...,m and consequently each plane of 

J r J 4 »J5>J6 

the sheaf of planes through the point P 4 is a determinantal locus for 


the satellite-position points. 

(iv) By equating either of the next three factors to zero, take 



one has the coplanarity condition of the points 

Pi = p 2 E (*2,y2> z 2)> P 3 E C x 3»y3> z 3)» Qj 4 5 (x j 4 » Y j 4 * z i 4 ) * 

Therefore the plane of the station-position points Pj , P 2 , P 3 is a 
determinantal locus for the set of the satellite-position points. Sim- 
ilarly from 

jp ^3 

A^ =0 or Aj =0 

one concludes that the plane (P^,P 2 ,P 4 ) is a determinantal locus for 
the set of the satellite-position points and from 



that this is the case for the plane (Pj,Pj,P 4 ). Notice that the last 
two loci (i.e, the planes (Pi,P2*P 4 ) an d CPi>P3>P 4 )) fall in the above 
considered case (iii). These are the determinantal loci for the satel- 
lite-position points independently from the relative position of the 
stations. 

Remark. In addition to the above determinantal loci one could quote 
those which are degenerate cases of the above. These determinantal loci 
lower further the rank of the jacobian matrix, when it has already been 
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lowered from six at least to five by one of the above quoted determin- 
antal loci. For example, if all the satellite-position points are on 
a straight line, this is obviously a determinantal locus that is inclu- 
ded in (32) for then there exists always a plane through P 4 and the 
straight line carrier of the satellite-position points. Such a loci 
which lower the rank further, once it has already been lowered from 
six, are trivial and result in an obvious way, being of no interest in 
this context. 

(2) Determinantal Loci for the Station-Position Points 

One is interested here for the determinantal loci for the station- 
position points independently from the satellite-position points. 
Since there are only four stations the only loci are the straight line 
determined by any two station-position points and the plane determined 
by any three of them. In the following the plane determined by any 
three station-position points is considered first. Suppose that the 
four station-position points are coplanar. Consider the determinants 

4 , 4 , 4 , 4 , 

expanded according to the first row, i.e, 

4 = a l x j + Pl Y j + Yl z j + 5 1 

a = 


A j 2 = a 2 X j + B 2 Yj + y 2 Zj + 

4 = a 3 X j + 0 3 Y j + Y 3 Z j + 6 3 

4 = a 4 X j + e 4 Y j + Y 4 Z j + 6 4 


(33) 1 
(33) 2 

(33)3 

(33), 


where ai, Yi> <5i , (i=l, 2,3,4) .are functions of the station coordin- 
ates only, (see relations 27). 
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By putting 


(34) 


Pi= e (a t + 


2 1/2 

Yi) ■ 


(e = sign of A^) 


one has 



P 1 p 2 P 3 p 4 

where d is the distance of Qj h (X.. ,Y. , Z - ) from the plane of the four 
station-position points. Indeed 

4 

-i ,(1 = 1*2,3, 4) 

Pi 

is the expression of the distance of the point Qj = (Xj,Yj,Zj) from the 
plane of three of the four station-position points, and if all four are 


coplanar (35) follows. 

From (35) and relations (E12) and (E13) of Appendix E one obtains: 

0 = a| - A^+ A^- A^ = d (p i-p2 + P3 _ P4.) ^ 
p l +p 3 = p 2 +p 4» ( d ^ °^ • 

Now consider the typical row of the jacobian matrix (Table 3) and effect 
the following rank-equivalent operations, in the sequence they are given 

(i) Divide each row by the corresponding A 4 (*) 

(ii) Subtract from x 2 -column, p 2 /p 4 times the x 4 -column. 

(iii) Add from x^-column, P 3 /P 4 times the x 4 -column. 

(iv) Add from yj-column, P3/P4 times the y 4 -column. 

(v) Change the signs of the columns x^, y 3 , x 4 , y 4 , z 4 > 

Then the rank-equivalent jacobian matrix is: 


*One may divide by any of A^ (i = 1 , 2 , 3, 4) 



x 2 

x 3 

y 3 

x 4 

y4 

z 4 

P? 

(X 2- X 4>p7 
— 

P 3 

(x 4“ x 3^ 

P3 

( >v^) P 7 

^r 

X 

1 

X 

v r y 4 

z r z 4 


Table 3. Rank-equivalent Condensed Jacobian Matrix for Four 
Coplanar Station-points. 

One observes that each of the first three columns has the same en- 
tries throughout all the rows. Therefore the determinants of all 6 x 6 
submatrices vanish and consequently the plane determined by any three 
of the four station-position points is a determinantal locus for the 
fourth one. 

By considering the determinants Aj , (i=l ,2,3, 4) one concludes that 
if three station-position points are collinear one of them vanishes. 

Thus the straight line determined by any two station-position points is 
a determinantal locus for each of the other two separately. 
Recapitulating the above results one has 
1 ) D eterminantal loci for the satellite-position points 

la) The four station-position points and any five satellite- 
position points determine through relation (30) a determinantal 
locus for the rest of the satellite-position points. 

lb) The plane through P 2 = (x 2 , Y2 » z 2 ) and perpendicular to 
the x-axis is a determinantal locus for the set of the satellite- 
position points. 

lc) Each plane of the pencil of planes through P 3 = (* 3 ^ 3 , 23 ) 
and perpendicular to the (x,y) -plane is a determinantal locus for 
the set of the satellite-position points. 
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l d) Each plane of the sheaf of planes through P 4 = (x 4 ,y 4 ,z 4 ) 
is a determinantal locus for the set of the satellite-position 
points. 

le) Each of the planes (P i ^2^3^ » ( p l p 2 p 4)> anc * (P1P3P4) a 
determinantal locus for the set of the satellite-position points. 
Although (P JP2P4) and (Pi p 3 p 4 ) are included in (Id) it seems appropriate 
to quote them along with (P1P2P3) which is not included in 

any of the above loci. 

2) Determinantal loci for the station-position points 

2a) The plane of any three of the four station-position points 
is a determinantal locus for the fourth one. 

2b) The straight line through any two station-position points 
is a determinantal locus for each of the other two separately. 

Remarks 

i) No loci that lower the rank of the jacobian matrix further, 
once it has already been lowered from six, is considered above. These 
loci are trivial and are easily obtained as further degeneration of 
those already given. For example, if the satellite-position points lie 
on a straight line, then there is always a plane through that line and 
the station-position point P 4 = (x 4 ,y 4 ,z 4 ), and this is nothing but an 
obvious further degeneration of the above quoted case (Id). 

ii) The above results correspond to the set of constraints 
X 1 j/l j z 1 >y2» z 2> z 3 = constant. However, for any other set of admissible 
constraints only the determinantal loci (lb), (lc), (Id), and (le) will 
be changed. Having the symmetric condensed jacobian matrix of Table 1, 
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one can carry out the above procedure for any other set of admissible 
constraints very easily. 

1.4.3 The Case of Five Ground Stations 

Let Pj s (*i»yi.*i). p 2 = ( x 2’ y 2 ,z 2^ P 3 F t x 3’ y 3’ z 3^ 

P 4 = (x 4 .y 4 .z 4 ), P 5 = (x 5 ,y 5 > z 5 ) be the position points of five observ- 

ing ground stations. Since the minimum number of participating stations 
at each event is four, there are 
N * C 4 + C 5 =5+ 1 = 6 

possible basic observational patterns, (see section 1.1.4). As a matter 
of fact, they may observe either all together or in groups of four. Let 
Qj = (Xj , Yj , Zj ) , (j = l , 2 , , . . ,m) be a set of satellite-position points 
observed in some pattern by the above five stations. Then one has a 


set of simultaneous functional relationships of the form* 

fij H (X r xi) 2 + (Y-j-yi) 2 + CZ^-Zi ^ 2 = r ly (37) 

This system contains 3m+3x5 = 3(m+5) unknown variables, six of which 
are to be constrained, so that there actually are 

p = 3(m+3) unknown variables. (38) 

Thus the number of simultaneous equations (37) must be 

v 2 u = 3 (m+3) (39) 

and for a unique least squares solution 

rank = p = 3(m+3) (40) 

where J M is the jacobian matrix of fjj with respect to Xj,Yj,Zj, 

x i ,y i ,z i’ (j = 1 > 2 »• ■ • i = l > 2 ,3,4,5) . If 

rank Jjq <p = 3(m+3) (41) 


*In this case for each j the index i takes on at least four out of the 
five values 1 , 2 , 3, 4, 5. 
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one has the determinantal loci of the configuration. In order to find 
out the determinantal loci one must search for the zeroes of the deter- 
minants of all yxy submatrices of the jacobian matrix J^. 

The expression of the jacobian matrix J M varies with the observa- 
tional pattern; however, it is not difficult to write down a general 
expression for the jacobian matrix, in the sense that it comprises the 
jacobian matrices of all possible observational patterns. In order to 
form that general expression, consider all possible station-quadruples, 
i . e. , 

0 5 E (1,2, 3,4) 

0 4 E (1,2, 3, 5) 

0 3 = d,2, 4, 5) (42) 

02 = (1, 3,4,5) 

0l i (2, 3,4, 5). 

Each satellite-position point is observed either by all five stations 
or by one of the above station-quadruples* It turns out that one can 
condense the jacobian matrix as in the case of four stations. The only 
new point in the procedure of getting the condensed jacobian matrix in 
this case, is that about the satellite-position points observed by all 
five stations. One verifies very easily by pivotal condensation that 
the condensed jacobian matrix corresponding to the case when all five 
stations observe simultaneously a set of satellite-position points, is 
formed by the condensed jacobian matrices of any two s tat ion -quadruples , 
say O5 = (1,2, 3, 4) and 0^ = (1,2, 3, 5) placed rowwise one after the 
other. It is understood, of course, that the condensed jacobian matri- 
ces of these two station-quadruples correspond to the very same set of 
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satel lite=position points. Hence it became obvious the way of forming 
the general expression for the jacobian matrix. Let {j,.}, {j^)> 

{j } {j } be the index-sets of the satellite-position points observed 
by the station-quadruples 05 , 04 , 0 ^, 02 , 0 ^ respectively, two and only two 
of which (no matter which ones) comprise* the satellite-position points 
observed by all five stations. By getting the condensed jacobian ma- 
trices of all station-quadruples separately, and placing them rowwise 
one after the other (regardless order) one forms the general expression 
for the condensed jacobian matrix for five stations. That is given in 
Table 4 where the symbol 

A ik’ i,k “ 1 , 2 , 3, 4, 5, 

is similar to A? that was used previously in the case of four stations, 
(see relations (27). Namely it stands for the determinant of 4-th 
order whose top row is occupied by the affine coordinates of the satel- 
lite-position point Qj (indicated by the upper index) and the rest three 
by the affine coordinates of three out of the five station-position 


points (in order of increasing index), excluding the stations Pj_ and 
(indicated by the lower indices). For example 



X- 

Y- 


J 

J 

43 ■ 

X 1 

y 1 

x 4 

f4 


x 5 

^5 

The order of the 

rows is 


zj 1 

2 4 1 

z 5 1 

conventional, though. 


(43) 


Some other order could 


be stipulated. 

*F.ach of those two station-quadruples comprises the satellite-position 
points observed by all five stations, in addition to those which are 
observed by that station-quadruple alone. 
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Table 4. The General Expression for the Condensed Jacobian Matrix for Five Stations in Range Networks 

























Having formed the general expression for the condensed jacobian 
matrix, one proceeds to find out the determinantal loci. At first an 
admissible set of constraints is chosen. As in the case of four sta- 
tions, take 

x l>yi> z l>y2> z 2* z 3 ~ constant. 

Since the columns of the jacobian matrix corresponding to the con- 
strained coordinates are zero, the condensed jacobian matrix corres- 
ponding to the above chosen set of constraints results from Table 4 by 
simply leaving "out" the columns corresponding to those constrained 
coordinates. That matrix has 15-6 = 9 columns. Thus in order to find 
out the determinantal loci one has to search for the zeroes of the de- 
terminants of all 9x9 submatrices. Again, as in the case of four sta- 
tions, the determinantal loci are distinguished in two categories: 
i) loci for the satellite-position points, and 
ii) loci for the station-position points. 

The determinantal loci of each category are independent from one 
another, in the sense that each one may occur independently from the 
others . 

1) Determinantal loci for the satellite-position points 
For a systematic exposition the following two groups of patterns, 
(comprising all the patterns) will be considered separately. 

(a) Two station-quadruples, say 0 5 = (1,2, 3, 4) and 0 5 e (1,2, 3,5) 
carry out all the observations. The theoretical treatment of this pat- 
tern includes that of the pattern when all five stations observe simul- 
taneously all the satellite-position points. 
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(b) Three, four or all five station-quadruples carry out all the 


observations . 

Remark. When all five stations observe simultaneously a certain set 
{Qj } of satellite-position points, one has the maximum number of equa- 
tions he can have with the same set {Q-j } and any possible observational 
pattern. The jacobian matrix corresponding to any observational pat- 
tern can be obtained from that of the pattern when all stations observe 
simultaneously by omitting certain equations. For example, suppose 
that the point-set 

{Qj}, je{j }={1,2,3, . . . ,m} 

has been observed by the quadruples 0^, O4 and O3 , so that O5 has ob- 
served the subset 


0 4 the subset {Qj }, j 4 e{j 4 ) C {jh 



where 

{j 3 }u{j 4 }u{j 5 Mjh ij 3 } n {j 4 -Cj 4 > n (j 5 {j 5 }n{j 3 }=<(>. 

The jacobian matrix of this pattern, results from the jacobian ma- 
trix of the pattern when all five stations observe simultaneously all 
the points in {Qj}, by omitting the rows corresponding to the observa- 
tions of stations P 3 , P 4 and P 3 toward the point-sets {Q- }, {Q^ }, and 

j 3 j 4 

{Qj } respectively. Therefore, by considering the jacobian matrix as 
if all stations observe simultaneously, no matter what the actual obser- 
vational pattern is, and checking its rank, if it is less than nine then 
this is the case for the jacobian matrix of the actual pattern as well; 
however if it is equal to nine, the rank of the jacobian matrix of the 
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actual pattern may or may not be equal to nine, and one must examine 
the rank of the actual jacobian matrix. 

From this remark becomes obvious why the pattern of two station- 
quadruples carrying out all the observations is so important to be 
treated separately, in searching for the determinantal loci in the case 
of five stations. 

(a) Two station-quadruples carry out all the observations 

Let 0 5 = (1,2, 3, 4) and 0 4 = (1,2, 3, 5) be the two station-quad- 
ruples. The condensed jacobian matrix of this case is given by the 
first two row-blocks of Table 4, 

{ V - {(X j 5 ’ Y j 5 ,Z j 5 )}> i 5 E{j 5 } ' 

is the set of the satellite-position points observed by 0 5 5 (1,2,3 ,4) 
and 

{Qi } 2 UX- , Y- ,Z, )}, j 4 e{j 4 }, 

J 4 J 4 J 4 J 4 

is the set of the satellite-position points observed by 0 4 = (1,2, 3, 5). 
Note that for 

(Q. } = {Q. ) = (Q 4 ) • (*) C 44 ) 

J 5 J 4 J 

one has the case when all five stations observe simultaneously all the 
satellite-position points, while for 

(Q.j >n'{ Q- } = <$>, («) 

3 5 4 4 

the case when there is no event at which stations P 4 and P5 enter to- 
gether. In the case when some satellite-position points are observed 
by all five stations, while the rest by either of the two station- 

*It is assumed that {Q, ) = {(L } is the set of all the satellite- 

J 5 J 4 J 

position points available. 
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Table 5. The Condensed Jacobian Matrix of the Case When the Station-Quadruples 0^ and O 4 carry out 
all the observations, under the Set of Admissible Constraints x-^, y-^, z^, y ^ z 2 > - 



























quadruples alone, neither (44) nor (45) is true. 

Considering the above adopted set of admissible constraints, the 
condensed jacobian matrix of the present case is given by Table 5 
(on preceding page), that is nothing else but the first two row blocks 

of Table 4 after omitting the columns of the constrained coordinates 
x p y^> y 2* z 2* z 3* 

For the determinantal loci one has to search for the zeroes of the 
determinants of all 9x9 submatrices - In order these determinants not 
to be identically equal to zero, each one of them must include at least 
three rows from each station-quadruple or that is equivalent at most 
six rows from each station-quadruple. 

Since the rank of the jacobian matrix is invariant under affine 
trnas format ion (see Appendix C) and the same is true for the determin- 
ants of its submatrices, one can express each of the determinants of 
the 9x9 submatrices as polynomials in the fundamental invariants (see 
Appendix F). This is done by Laplace expansion in the same way that 
was done for the case of four stations. Thus the determinant, 
|a( 0 5 ,0 4 )|, of the submatrix corresponding to a certain selection 
kpk ? , . . . ,k g of satellite-point-indices (satisfying the following con- 
ditions (47) is given by the following messy, but nevertheless useful, 
summation formula. 
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|a(o 5 ,o 4 )| 

P 1 (p2( p 3^) 



x 3 /3 1 

*2 1 


h 2 \ 2 ' 

X kl 1 


X k 3 , kj 1 


k, 

* 2 \ 


k-> k-» 

fl 3u 6 3 v 


where 


k 4. k 5, k 6,, k 7, k 8, k 9 k 4 k 5 k 6 k 7 k 8 k 9 
7*4 *>5 fi 4S 6 4S & 4S A S4 4 b4 A 54 


(46) 


ki e { j 5 } U {^4} » 
k 2> k 2 e { 35 } U { j 4 > , 
k 4, k 5 , k 6 e » 

k 7> k 8» k 9 e ^ 4 ^ > 


(47) 


\,u,v 


4, if kj, k 2 , k 3 z (jg) , respectively, 

5, if k^, k 2 , kg £ { 34 } , respectively. 


(48) 


^k ,k ,k 
X»4 


x 4 y 4 z 4 1 



, p k 7> k 8» k 9 = 


x 5 y& z 5 1 

Xk 7 Yk 7 Zk 7 1 

x k Y i- z v 1 
k g k g kg 

X. Y, Z. 1 
k 9 k 9 k 9 


(49) 


By equating the determinant of a 9x9 submatrix to zero, i.e., 

[a(0 5 ,0 4 ) I = 0 (50) 

one has a hypersurface with at least 24 and at most 36 variables. One may 
keep constant the coordinates of all the points involved but one, and consider 
(50) as a three dimensional locus of that point. However the form of (50) 
varies with the choice of the submatrix in such a way that searching for the 
three dimensional second order determinantal loci is not as simple a task as 
in the case of four stations. Before proceeding systematically, take for 
example the determinant of the following Table 6. 
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Table 6 - A 9x9 Submatrix of the Table 5. 



The above table is a 9x9 submatrix of the jacobian matrix, which has 
five rows from the station-quadruple 04=(1,2,3,4) and four rows from 
Oc (1, 2, 3,3). By considering its determinant as a functional with re- 
spect to X 8 , Yg, (the coordinates of Q^-one of the points observed 
by 0$) and equating it to zero, one obtains a second degree surface as 
the locus of the point Qs= (X5 .Y^, Z5) , that passes obviously through 
Ql>Q2'Q3>Q4> Pl»^2»P3>P4 anc * t * i e origin (0,0,0). For each selection 
of the last four rows in the station-quadruple 04=(1, 2, 3, 5) the deter- 
minant of the above table equated to zero, gives a second degree sur- 
face as the locus of Q5=(X5,Y5, Z5) that passes through the same nine 
points. Provided these nine points do not constitute a singular con- 
figuration, they determine uniquely a second order surface passing 
through each one of them. Therefore the loci of Q5- (X5 ,Yg ,Z 8 ) (second 
degree surfaces) corresponding to all possible selections of four rows 
in 04^(1, 2, 3, 5) coincide with the second degree surface determined by 
the nine points Qj,Q 2>Q3>Q4> p i> p 2> p 3>^4 an< * t * ie ori g in > provided these 
points do not constitute a singular configuration. Note that instead 
of Q g one may consider any of the satellite-position points observed by 
0 5 =(1» 2, 3,4) . 

Similarly, by considering the determinant of the above Table 6 as a 
functional with respect to X g ,Y 9 ,Z 9 (the coordinates of Qg-one of the 
points observed by 045(1,2,3,5)) and equating it to zero, one obtains a 
second degree surface as the locus of Qg=(Xg ,Yg,Z 9 ) , that obviously 
passes through the points Q 6 ,Q 7 ,Q 8 , P1.P2tP3.P5 and the origin (0,0,0). 
For each selection of the first five rows in the station-quadruple 
°5 e ( 1> 2 »3,4) the determinant of the above table being considered as a 
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functional of X 9 ,Y 9 ,Z 9 and equated to zero, gives a second degree sur- 
face as the locus of QgHCXg.Yg, Z g ) that passes through, the same eight 
points. But there are infinitely many second degree surfaces through 
eight points. Therefore the loci of Q 9 =(X 9 ,Y 9 ,Z 9 ) (second degree sur- 
faces) corresponding to different selections of five rows in 
0^=(1,2,3,4) are generally different. However, instead of Q 9 = (X 9 ,Y 9 ,Z^| 
one may consider any other of the satellite-position points observed by 
0 4 e( 1, 2,3,5) . That is, the second degree surface that is , determined by 
any selection of five rows in (1 , 2, 3,4) , say those corresponding to 
Ql ,Q 2 ,Q 3 ,Q 4 >Q 5 » and an y selection of three rows in 0 4 =(1,2,3,5), say 
those corresponding to Q6*Q7>Qs» passes through all the satellite- 
position points observed by the station-quadruple ,0 4 = (1 , 2, 3, 4) , through 
the station-position points Pi»P 2»^3’^5’ and or ^8^ n * Therefore, if 
the station -quadruple 0 4 =(1, 2, 3,5) , observes more than three satellite- 
position points, the very same second degree surface that is deter- 
mined by the nine points: origin, Pj , p 2 » P 3» P 5* and an X f° ur satel “ 

lite-position points from the station -quadruple 0 4 =(1,2,3,5) , it is 
determined completely also by Pj,P 2 ,P 3 ,Ps, origin, and any three satel- 
lite-position points. Of course this is not the case in general with a 
second degree surface. 

In order to find out the determinantal loci of second degree, one 
may proceed as follows. , 

As it was pointed out before, each 9x9 submatrix being considered 
for the determinantal loci must have at least three rows from each 
station-quadruple and at most six rows; otherwise its determinant is 
identically zero. Thus one considers all 9x9 submatrices of the 
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jacobian matrix falling in four categories: 

(а) The submatrices which have six rows from the s tat ion -quadruple 
O5H (1, 2,3,4) and three rows from O4S (1 , 2, 3, 5) , 

(0) The submatrices which have five rows from 052(1,2,3,4) and four 
rows from 0^2(1, 2,3, 5), 

(X) The submatrices which have four rows from 0^2(1, 2, 3, 4) and five 
rows from 0^2(1 ,2, 3, 5) , and 

(б) The submatrices which have three rows from 052(1,2,3,4) and six 
rows from 0^2 (1,2, 3,5) . 

The determinantal loci for the satellite-position points observed by 
the station-quadruple 0 S 2(1,2,3,4) will be considered below. 

The determinantal loci for the satellite-position points observed 
by the station-quadruple O42 (1,2, 3, 5) result in exactly the same way. 
One simply interchanges the roles of 0 5 and 0 4 in the following discus- 
sion for the determinantal loci of the points observed by 0 S 2(1, 2, 3,4) . 

for the rank of the jacobian matrix to be less than nine the deter- 
minants of the submatrices of all the categories ((a), (6), (Y) , (6)) 
must vanish simultaneously. Thus the loci for the satellite-position 
points observed by 052(1,2,3,4) resulting from the four categories must 
coincide, since they refer to the same set of points. 

Prom the category (a) : Let Qj ^2 ^3^4 ^5^5 be six points from 

the station-quadruple 052(1,2,3,4) and Q 7 ,Q g ,Q 9 from 0 4 =(1,2,3,5). By 
considering the determinant of the submatrix corresponding to these 
points as a functional with respect to the coordinates of one of the 
points observed by 052(1,2,3,4), say Q^2 (X^,Y^, Z^) and equating it to 
zero, one has a second degree locus for »^6 ’^ 6 ^ ’ passes 
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through the points Qj , Q2>Q3»^4 ^5’ P 1* P 2* P 3» P 4 and the ori 8 in 
For each selection of three rows in 04 = ( 1 , 2 , 3,5) the locus of 
Q 6 h(X 6 ,Y 6 ,Z 6 ) passes through the same points.' Since this is a second 
degree surface, it can be determined uniquely by any nine of its points, 
provided they are not in singular configuration. Moreover instead of 
Q 6 =(X 6 ,Y 6 ,Z 6 ) one may take any other point from those observed by 
0 5 =(1,2, 3,4). Thus the locus of the satellite-position points observed 
by 0 5 =( 1 , 2 , 3 , 4 ) , dictated by the category (a), is a second degree sur- 
face whose equation can be given either by equating to zero the deter- 
minant of a submatrix of the category (ot) ! (being considered as a func- 
tional with respect to the coordinates of one of the points observed by 
0 5 =( 1 , 2 , 3 , 4 )) , or by determining the second order surface which passes 
through the nine points: origin, Pi>P 2 ,P 3 ,P 4 and an ^ ^ our satellite- 

position points, observed by 05 = (1 , 2, 3, 4) . Let this locus be denoted 
3 

by L 5 . 

Fro m the category (3) : This category was considered in the above 

given example. It was found that, for each selection of four rows in 
045 ( 1 , 2 , 3 , 5 ) the locus of the satellite-position points observed by 
0 S =(1,2,3,4) passes through the points: origin, Pi,P2* p 3» p 4 and four 

from the points observed by 055 ( 1 , 2 , 3 , 4 ). This locus is again a gener- 
al second degree surface whose equation is given either by equating to 
zero the determinant of a submatrix of the category (3) , (being consid- 
ered as a functional with respect to the coordinates of one of the 
points observed by 0^(1 , 2, 3, 4) ), or by determining the second order 
surface which passes through the nine points: origin, Pi»P2> p 3» p 4 and 

any four points observed by 055 ( 1 , 2 , 3 , 4 ). Obviously this locus 
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3 4 

coincides with L 5 of category ( a ) . Let it be denoted by L5. 

From category (Y) : This category was also considered in that, 

there (second part of the example) the discussion was for the locus of 
the points observed by 0^(1 ,2, 3, 5) and not 0 4 =(1 , 2, 3, 4) . However, the 
results are similar. Thus by transferring those results here one has 
that, for each selection of five rows in 0 4 =(1,2,3,5) the locus of the 
satellite-position points observed by 0g=(l,2,3,4) passes through the 
points: origin, P^,P ? ,P2,P 4 and any three points observed by 

0^=(1,2,3,4) . But eight points do not determine a general second degree 
surface, i.e., there are infinitely many general second degree surfaces 
passing through eight points. However, all the loci corresponding to 
different selections of five rows from 0 4 e( 1,2,3,5) refer to the very 
same set of points that is the points observed by 0cj=(l,2,3,4) . There- 
fore they must coincide with Lg and L^, in case of a determinantal 
locus. That implies that the locus of the satellite-position points 
observed by 0^=(1, 2,3,4) is not a general second degree surface, but it 
is such that any eight of its points determine it completely. The equa- 
tion of this locus is given by equating to zero the determinant of a 
submatrix of the category (y), (being considered as a functional with 
respect to the coordinates of one of the points observed by 0g=(l ,2, 3,4) ), 
For the existence of such a locus take the equation of b| (=l£) and 
verify that the same locus L s is determined completely by any eight of 

5 

its points. Let such a locus be denoted by L 8 . 

From the category (6) : Let Qi,Q 2*Q3 be three points observed by 

0 5 e( 1,2,3,4) and Q 4 ,Q 5 ,Q 6 ,Q 7 ,Q 8 ,Q 9 six points observed by 0 4 = (1 ,2, 3 ,5) . 
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By considering the determinant of the submatrix corresponding to these 
points as a functional with respect to the coordinates of one of the 
points observed by Ot^ (1, 2, 3, 4) , say , and equating it to 

zero, one obtains a second degree surface as the locus of Qj^CXjjYjjZ^) , 
that passes through the points Qi»Q 2»^ 1»^2 , ^3 , ^ > 4 an< ^ or ^S^ n =(0,0,0). 
For each selection of six rows in 0^(1, 2, 3, 5) the locus of 
Qj=(Xj, Y j, Zj) passes through the same seven points. Instead of Qj one 
may consider any point observed by 0tj=(l,2,3,4) . Thus the locus of the 
satellite-position points observed by 05=0,2, 3,4) in this category, 
passes through the points: origin, PijP 2'^3*^4 anc * an Y two ^ rom t ^ ie 

points observed by 0^= (1, 2, 3, 4) . The equation of this locus is taken 
by equating to zero the determinant of a submatrix of this category (6), 
(being considered as a functional with respect to the coordinates of 
one of the points observed by 0gS(l,2,3,4)) . This locus must coincide 
with the loci L5, L*, and L5 of the previous categories, in case of a 
determinantal locus. That implies that the locus of the points observed 
by C>5=(1,2,3,4) is not a general second degree surface but it is such 
that any* seven of its points determine it completely. Let such a locus 
be denoted by L^. For the existence of such a locus take the equation 
of L5 (=1*5) and verify that the same locus is determined completely 
by any six of its points. 

The same discussion may be repeated for the determinantal loci of 
the satellite-position points observed by 0^(1, 2, 3,5) , to find the com- 
pletely analogous loci L4, L4, L^, L^. In the general case, when both 
station-quadruples 0g=(l,2,3,4) and 0^=(1,2,3,5) observe six or more 
points each, the determinantal locus is L5 for the points observed by 
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0 5 =(1,2,3,4) and L 4 for the points observed by 0 4 e ( 1 , 2, 3, 5) . It is 
understood that these loci, and L 4 , go together. That means that 
for the rank of the jacobian matrix to be less than nine, both the 
satellite-position points observed by C> 5 =( 1 , 2 , 3 , 4 ) must lie on Lc-, and 
the satellite-position points observed by 0 4 e( 1,2,3,5) must lie on L 4 . 

It is possible l| and l| to coincide. That happens whenever they have 
seven* points in common; and since they do have in common the points 
* > 1 , * > 2 , * > 3’* 3 and the or igi n » at least three satellite-position points 
observed by all five stations are sufficient for the loci L 5 and L 4 to 
coincide. Besides the general case when each quadruple observes six or 
more points, one has some specific cases, which are those when one of 
the quadruples observes either only three or only four or only five 
points. Obviously, if one of the quadruples observes only three points, 
one has to consider only the category (a), if it observes only four 
points, one has to consider only the categories, (a) and (B), and if it 
observes only five points, one has to consider the categories (a) , ( 0 ) 
and (Y) . 

After the second degree determinantal loci, one proceeds to con- 
sider the linear ones. They may be taken from the expression (46) by 
equating to zero each factor of the typical term of the summation. Thus 
one has: 


(i) From 


x 2 

*kj 


= 0, (kie{ J5} U { j^}) 


"Seven points, ^because as it was said above (category ( 6 )) L 5 (and 
counterlocus L 4 ) is not a general second degree surface, but it is 
that any seven of its points determine it completely. 


(51) 

its 

such 
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it is concluded that if all the satellite-position points, (observed 
by both station quadruples), lie in a plane perpendicular to the x-axis 
and passing through the point P 2 S C X 2 »y 2 » z 2 ^ > *he column of all 

9x9 submatrices i a zero column. Therefore that plane is a determin- 
antal locus for all the points. 

(ii) From 

= 0, (k 2 ,k 3 e{ j 5 > U {j 4 }) C 52 ) 

it is concluded that if all the satellite-position points lie in a 
plane through the point P 3- ( x 3 > 73 * * 3 ) ar *d perpendicular to the (xy)- 
plane, the determinants of all 9x9 submatrices are zero. Therefore 
each plane of the pencil of planes with axis the projection line of P 3 
to the (x.y)-plane is a determinantal locus for all the points. 


x 3 73 

S y k 2 


**3 Y k3 1 


(iii) From 


k 4 ,k 5 ,k 6 

J>4 


x 4 y 4 z 4 

X k 4 Y k 4 Z k 4 
x k 5 *k 5 z k5 
x k 6 Y k 6 z k 6 


1 

1 

1 

1 


(k 4 ,k 5 j k 6 e{j 5 }) 


(53) 


it is concluded that, if the satellite-position points observed by the 
station-quadruple 0^(1, 2, 3,4) lie in a plane through the point 
P 4 = (X 4 , 74 , z 4 ) j then the determinants of all 9x9 submatrices are zero. 
Therefore, each plane of the sheaf of planes through the point P 4 is a 
determinantal locus for the satellite-position points observed by 



(iv) From 


x s ys. z 5 1 

x k7 Y k7 Z k7 1 

x k 8 % z k 8 1 

% % s i 

it is concluded that, if the satellite-position points, observed by the 
station-quadruple 0 4 e( 1, 2,3,5) , lie in a plane through the point 
p 5=( x 5>y5> 25) > then the determinants of all 9x9 submatrices are zero. 
Therefore each plane of the sheaf of planes through the point P^ is a 
determinantal locus for the satellite-position points observed by the 
station -quadruple 0 4 =(1, 2,3,5) . 

(v) From either one of the following three relations 

&45 = 0, =0, = 0,(k4,k5,k6e{j5» , (5S) 

it is concluded that the plane ( p i, p 2> p 3) is a determinantal locus for 
the set of the satellite-position points observed by 053(1,2,3,4). 

(vi) From either one of the following three relations 

a 54 ~ °» A 54 “ a 54 = °» (k7,k8,k9e{j4)), C 56 ) 

it is concluded that the plane ( p i, p 2» p 3) is also a determinantal locus 
for the satellite-position points observed by 0 4 =(1, 2,3,5) . 

Note: In the last two cases, (v) and (vi) , it is observed that, one may 

have the same determinantal locus for different sets of satellite- 
position points. 

(vii) From 

ki . 4 > if ♦k^e{j 5} 

A 2A = °» U { 5 4 } » A = (57) 

J 5, if k 1 e{j 4 ) 


_k7 ,k8,k9 
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it is concluded that, if the satellite-position points observed by 
O s =(l,2,3,4) lie in the plane (P^^s) and simultaneously the satel- 
lite-position points observed by 0 4 ^ (1 , 2, 3,5) lie in the CPi» p 3 »P 4 )» 
then the determinants of all 9x9 submatrices are zero. Therefore the 
planes and (P 1 ,P 3 ,P 4 ) are simultaneously determinantal loci 

for the satellite-position points observed by 0 5 =(1,2,3,4) and 
0 4 = (1 , 2,3,5) respectively. 

(viii) From either one of the following two relations 

k 2 k3 k2,k3e{j s)U{ j 4 ), C 58 ) 

a 3u = °> a 3v = . _ 4, if k 2 ,k 3 e{j 5 > respectively 

" 5, if k2, 3 ed 6 ) respectively 

it is concluded that the planes (Pi,P2>P4) an ^ (Pl»P2>^5) together con- 
stitute a determinantal locus for the point-sets {J 5 } and {j 3) respec- 
tively. These are the determinantal loci of the satellite-position 
points observed by two station-quadruples. 

(b) Three, four or five station-quadruples carry out all the observa- 
tions . 

It was said above that, if a determinantal locus of the previous 
case of two station-quadruples occurs, then that is also a determinantal 
locus for any other observational pattern with the same satellite-posi- 
tion points. But if the previous case is non-singular, (i.e., its con- 
densed jacobian matrix is of rank nine) , that does not mean that this is 
the case for the other patterns with the same set of satellite-position 
points. Thus one has to search for the determinantal loci of the other 
observational patterns too. 

These cases are treated in a similar way. For the determinantal 
loci of second degree one may partition the set of all 9x9 submatrices 
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of each of these cases in categories as before, while the linear deter- 
minantal loci are taken again from the expressions of the determinants 
of all 9 x 9 submatrices as polynomials in the fundamental invariants, (see 
Appendix F) . 

Consider the case where three station-quadruples carry out all the 
observations. Without loss of generality take the quadruples 
0 5 e(1,2,3,4), O 4 E (1,2, 3,4) and O 3 E ( 1 ,2, 4, S) . The jacobian matrix of this 
case is given by the first three row blocks of Table 4. For the set of 
admissible constraints 

x l l)r l' z l ,y 2 ’ z 2 ’ z 3 = constant, 

leave out the columns corresponding to these coordinates. For the deter - 
minantal loci of this case one has to examine the zeroes of the determin- 
ants of all 9x9 subma trices. Each 9x9 submatrix of this case has at 
least one row from each station -quadruple, and it must have at least 
three rows from any two out of the three quadruples, in order its deter- 
minant not be identically equal to zero. Consider as an example the 9x9 
submatirx on the following page. 

bet {j 5 }, (34}. (j 3 } be the sets of the satellite-position points 
observed by the station-quadruples 05 =( 1 , 2 , 3 , 4 ) , 04 =( 1 , 2 , 3 , 5 ), 

04 = 0 , 2 , 4 , 5) respectively. 

Second degree determinantal loci . Working as in the case of two 
station-quadruples one finds for the general case, (when each one of 
the three station-quadruples observes six or more satellite-position 
points), that the second degree determinantal locus consists of three 
second degree surfaces L^.L^L^, one for each of the sets {I 5 }, (34 3, 

(j 3 } respectively. The locus L 5 for the points ( j 5 }, passes through 
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Table 7. A 9x9 Submatrix of the Jacobian Matrix of the case when the three station-quadruples 
0 5 =(1,2,3,4), 0 4 =(1,2,3,5) and 0 3 e( 1, 2,4,5) carry out all the Observations. 



the points: origin, P 1 .P 2 .P 3 , and p 4 , and its equation is taken by 

equating to zero the determinant of any 9x9 submatrix which has only 
one row from 0 5 = (1 , 2 , 3 , 4) , say the corresponding to Qj s E ( x j 5 » Y j 5 » z j 5 ) > 
that determinant being considered as a functional with respect to 


Xj Yi Z- . This is a second degree surface which contains the points 
J5’ J5’ J 5 

{j 5 >. Similarly the locus L 4 for the points { j 4 } , passes through the 
points: origin, Pi,P 2 . p 3 » and P 5 , and its equation is taken by equating 
to zero the determinant of any 9x9 submatrix which has only one row from 


0 4 =( 1 , 2 , 3 , 5 ), say the corresponding to Qj 4 = ( x j 4 , Y j 4 ,Zj 4 ) , that determin- 
ant being considered as a functional with respect to x j 4 » Y j /( » Z j 4 ’ At 
last the locus L 3 for the points { j 3 } , passes through the points: 
origin, p i, p 2 » p 4 > and p 5 » and its equation is taken be equating to zero 
the determinant of any 9x9 submatrix which has only one row from 


0 3 =(1, 2,4, 5) , say the corresponding to Q j C^j ' that determin “ 
ant being considered as a functional with respect to 3 ,Zj 3 . Argu- 

ing as in the case of two station-quadruples, one concludes that each of 
the loci L 5 ,L 4 ,L 3 is not a general second degree surface, but it is such 
that any five of its points determine it completely. Since in pairs 
these loci have four points in common, if there exists a point observed 


by all five stations, those loci have to coincide with each other. 


Linear determinantal loci . The expression of the determinant of any 
9x9 submatrix of the jacobian matrix of this case, as a polynomial in 


the fundamental invariants is the following mess. 
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where 

Hjef j 5 > U( J4} U( J3}, 

n 2 ,n 3 e{ j 5 }|j{ 14^. 

n 4,n 5 ,n 6 e{ j5>U { 33), 

n 7> n 8’ n 9 e {^4}U{ j e K 

< 5 , if n 4 ,n 5 ,n 6 e( J5} respectively, 
(3, if n 4 ,n^e{ j 3 ) respectively. 


: 1 » K 2> 


x l> x 2> x 3 = 


M = 


4, if ny,ng,nge{ J 4 } respectively, 

3 , if n 7 ,ng,rig£{ J 3 } respectively, 

5, if njcfjs}, 

4, if njc{ J 4 }, 

3, if nje{j 3 }. 


( 59 ) 


V 1 » v 2 “ 


(5, if n 2 ,n 3 e{ J 5 } respectively, 
i4, if n 2 ,n 3 e{j 4 } respectively. 


and ^)^ n 4» n 5> n 6^^ n 7> n 8> n 9 are gi ven by ( 49 ). 

By equating to zero each factor of the typical term in the expression 
(59) one obtains the linear determinantal loci. 
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(i) From 

= 0, (njdj 5 > U {j 4 > U f j 3 >) ( 6 °) 

x m i 

one concludes that the plane through the point = ( x 2 »>'2 > z 2' and perpen- 
dicular to the x-axis is a determinantal locus for all the points 

{j 5} iHj 4^ U(j 3} • 

(ii) From 

x 3 x 3 1 

x n2 Y n2 1 =0 (n 2 .n 3 £{j 5} U {j 4 > (61 ) 

X n3 Y n3 1 

one concludes that each plane of the pencil of planes with axis the pro- 
jection line of P3=(x 3 ,y3,z 3 ) to the (x,y)-plane is a determinantal 
locus for the point-set 


(iii) From 

< 3 > 4 4 ’ n5 ,n& = 0 , (n 4 ,n 5 ,n 6 e{j 5 } 3 ^) ( 62 ) 

one concludes that each plane of the sheaf of planes through the point 

P 4 =(x 4 ,y 4 ,z 4 ) is a determinantal locus for the point-set { j 5 > U tj 3 > - 

(iv) From 

pn 7 ,n 8 ,,n 9 = ^ (n7>n8>n9 { j 4 }(j{ j 3 }) (63) 

one concludes that each plane of the sheaf of planes through the point 

P5 = (x 3 ,y 3 , z 3 ) is a determinantal locus for the point-set (j 4 }U{j 3}- 

(v) From either one of the relations 

n 4 ,n s ,n 6 e(j slUti 3} ( 64 ) 



5, if n 4 ,n 5 ,n 6 e{j 5 irespectively 

K 1 * K 2 * K 3 = 

3 , if n 4 ,n5,n 6 e{j 3)respectively 
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one concludes that the two planes (Pi, p 2> p 3) ancl ( p l> p 2> p s) constitute 
simultaneously a determinantal locus for the point-sets (35) and { j 3 } 
respectively. 

(vi) From either one of the relations 

n 7> n 8» n 9 c ^ 4} U {j 3} (65) 

1 4, if n^ng.ngeU 4)respectively, 
3, if n 7 ,n 8 ,n 9 e{j 3>respectively 
one concludes that the two planes ( p l, p 2> p 30 and ( p i, p 2» p 4) constitute 
simultaneously a determinantal locus for the point-sets U4} and {j 3} 
respectively. 

(vii) From 

njeij 5} U(j 4} U(j 3^ 

5, if njcCjs) 

4, if np:p4} (66) 

3, if nic{j 3} 

one concludes that the planes ( p l» p 3> p 4)> ( p l> p 3« p 5) an< ^ ( p l> p 4> p 5) are 
simultaneously a determinantal locus for the point-sets {35}, (j4 )j (j 3} 
respectively. 

(viii) From either one of the relations 

n 2 i n 3 e (i 5 } U ( j 3} (67) 

1 5, if n2,n3e{j 5} respectively 
4, if n2,n3e{j 4} respectively 
one concludes that the planes (Pi, p 2, p 4) and (Pi,P2, p S) constitute 
simultaneously a determinantal locus for the point-sets {35} and {34} 
respectively. 




n 7 ns ng 

A 5Xi = °> ^5A 2 = °» A 5A 3 = 0 
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For the cases when four or five station-quadruples carry out all 
the observations, one proceeds in a similar way. 

2) Determinantal loci for the station-position points 

The only loci determined by the station-position points alone, in- 
dependently from the satellite-position points are: 

(a) The straight lines determined by pairs of station-position 
points, 

(b) The planes determined by triples of station-position points, 

and 

(c) The sphere determined by any quadruple of station-position 
points . 

The sphere determined by any quadruple of station-position points 
is not a determinantal locus independent from the satellite-position 
points, as it was seen above while treating the determinantal loci for 
the satellite-position points. One can realize what would happen if the 
sphere of any four station-position points were a determinantal locus for 
the fifth one, independently from the satellite-position pointsl Then 
the configuration would be near singularity by fact; since the surface of 
the earth on which the stations lie is approximately spherical. 

For the determinantal loci consisting of straight lines, take first 
the observational pattern when two station-quadruples, say 05=(1,2,3,4) 
and 04= (1 , 2, 3 , 5 ) , carry out all the observations. 

Consider the determinants 



in the condensed jacobian matrix of Table 5 . 
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From the determinants 

a 2S* a 24 C for an y i'* 

of the X 2 -column block it is concluded that, if each of the point-ti iples 
t p i- p 3- p 4 } .»' {P 1’ P 3> P 5 } 

is collinear, the x 2 -column becomes a zero one. Therefore the straight 
line (Pj^Pj) is a determinantal locus for P 4 and P 5 . 

Similarly from the rest column blocks of the condensed jacobian 
matrix of Table 5 one obtains the following. From the x 3 ,y 3 -column block 
it is concluded that* if each of the point-triples 
{p 1j p 2 ,p 4 }, {Pi.P2.P5} 

is collinear the x 3 ,y 3 -column block becomes a zero one. Therefore the 
straight line (P 1 .P 2 ) is a determinantal locus for P 4 and P 5 . 

From the x 4 ,y 4 , z 4 -column block and X 5 >y 5 » z 5 _c °l umn block, it is con- 
cluded that, if the point-triple 
{Pl.P2. P 3> 

is collinear both of the last two column blocks become zero ones. There- 
fore the straight line determined by any two of the three points Pi,P2> p 3 
is a determinantal locus for the third one. 

These are the determinantal loci consisting of straight lines for 
the specific observational pattern when the two station-quadruples 
0 5 =(1 2,3,4) and 0 4 =(1,2,3,5) carry out all the observations, and under 
the set of admissible constraints x^ ,y\, z i,y2» z 2> z 3 = constant. Note 
that, if at least one of the coordinates x^,yj,zj of P]^ were unconstrained, 
one would have also the straight line (P 2 .P 3 ) as a determinantal locus of 
P 4 and P 5 . Thus the selection of the set of admissible constraints does 

affect this type of determinantal lock. In the following a general pro- 
cedure of determining such a determinantal locus is described. This 
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procedure is dictated by the form of the condensed jacobian matrix and 
it is applicable to any number of stations and any observational pattern. 

It is observed that each determinantal locus of this kind is associated 
with a station. In other words there are as many such loci as stations. 
Consider any station of the network, say • The determinantal locus 
associated with this station is determined as follows. Let Oj , 0 2 , 0 ^ 

be the station-quadruples where Pj_ participates in. Bv excluding 
there remains from each one of the station-quadruples 0j,...,0k a point- 
triple. Let these triples be denoted by 

Pi.P2.P3 - Pl.Pi.P| p l. p 2. p 3 • 

Let also 

(c 1 * (e 2) » • ■ • > (e k) 

be the straight lines determined by any two of the three points of the 
above point-triples respectively. Without loss of generality one may 
take 

fe^ECP^P^ , fe 2 ) = (Pi.P2)> * ’ fek) 5 ( p i> p 2) 

Then the straight lines (ej), (e 2 ), constitute simultaneously 

12 k 

a determinantal locus for the points P3,P3, P3 respectively. 

For the case of five stations, since any two station-quadruples have 
three station-position points in common, and consequently any two point- 
triples from the above have two points in common, one must have: 

fe 1 ) = (e 2 )= ••• 

For the determinantal loci consisting of planes, one observes that, 
if any two station-quadruples are coplanar, then their planes coincide, 
for they have in common three station-position points. Therefore one may 
proceed as follows. 
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Suppose that the five station-position points are coplanar. 
Consider the determinants 

A {k> (i,k= 1,2, 3,4,5) 
expanded according to the first row, i.e. 

A ik = «ik*j + 0ik Y j + Yik z j + 6 ik» ( 68) 

where aik, Bik, Yik> *ik are functions of the station-coordinates only , 
(corapare with relations (33). 

By putting 

p ik = e(a? k + 3- k + Y-^ 1/2 < £Esi S n of A ik 5 (69) 

one has (compare with (34) and (35)) 


1 

Aik 

Pik 


= d, (for all i,k and j) 


(70) 


where d is the distance of Qj _ (Xj , Yj ,Zj ) from the plane of the station- 
position points. 

Now consider the general expression for the condensed jacobian 
matrix, Table 4, along with the set of admissible constraints introduced 
previously, (i.e., , Y\ , Z\ ,y? » z 2 > z 3 = constant), and effect the follow- 

ing rank-equivalent operations. 

(i) Divide each row j 5 of 0 5 "(1,2,3,4) by the corresponding 

A j S 

*45 . 

(ii) Divide each row j 4 of 04 =( 1 , 2 , 3 , 5 ) by the corresponding 

A j4 

*54 . 

(iii) Divide each row of 0^(1, 2,4,5 by the corresponding 

j 3 

*23. 

(iv) Divide each row 12 of (1,3, 4, 5) by the corresponding 

4i 
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(v) Divide each row jj of 0-^(2, 3,4, 5) by the corresponding 


(vi) Subtract from X2-colurnn, P25/P45 times the X4~column. 

(vii) Subtract from X5-column, P24/P54 times the x^-column. 
viii) Add to x 3 -column, p 35 /,f> 45 times the X4 -column. 


(viii) Add to x 3 -column, 
(ix) Add to yj-column, 
(x) Add to xj-column. 


p 3 s / p 45 times the y 4 ~column. 
P 34 /P 54 times the xg-column. 
P34/P54 times the y^-column. 


(xi) Add to y 3 -column, P34/P54 times the y 5 -coiumn. 

(xii) Subtract the first row in each station-quadruple from the 
rest of the rows in it. 


Then it is obtained: 
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0 


0 


Xj 41 ” x j 42 | Y j 4 l" Y j 42 z j 41 * z j 42 


(XjJl-X4)^f ]W}Jl-y4)^| |( Z )31-IS^ ( X j31~*5^ j< Y j31-r5>rej - 

CXp2- x j3I^!t Y j32- Y )31 ^i (Z j 3r Z J3I^J^j32- x j3l)^i^j32-^J3 l)^|| < Z j32- z j3tVff 


C2 (Xj21'X4)^| i(Vj 21 -y 4 )^| ( z j 2 i-* 4 )^f CXjn-^gl |tVj2i-ys)g§ Ujzi-zsyff 


( x <22- x i2l)^^ Y J22- Y J2l)^Tlt Z i22- Z }21^|t x j22- X J21>4l (Y J22- Y j21>P32 (Z 3^~ Z i zl) g37 



<Vr x 4^ < z m- x 4^ i (x nr*5>gfr pir^Sfr ^n-VSfr , 
(x ji2‘ x jn ) ?if (I ii2- I in ) ^f! cx ii2' x jn 1 gfr tY ji2~ Y j 11 j g |f 


The general expression for the condensed jacobian matrix for five 
coplanar stations, transformed by rank-equivalent row-operations 
under the set of admissible constraints x^yp ,y 2 , z 2 > z 3 “ constant. 












Note . After effecting the first eleven operations, each of the first 


three columns in each station-quadruple, has the same entries through 
out all the rows in that quadruple. Thus after the twelfth operation, 
the first three columns have zeroes everywhere except in the first row 
in each station-quadruple. 

One verifies immediately that any three of the five non-zero rows 
of the first three columns in Table 8 are linearly dependent. Consider 
for example the first three non- zero rows, i.e.. 


*1 =[( x 2* x 4)p“> ( x 4“ x 3)p||-» 04-y^py||] 
^2 =i ( x 2 _x 5)|r§^ t X 5 _x 3)pr||, (/s-J^pi^-] 


^3 =[ a 3 . 


t>3> 


c 3 


( 71 )! 

(71)2 

(71)3 


where 


a- = xo + y , p 43 p 25 . _ p 53 p 24 
3 2 4p 23 p 45 Sp 23 p 54 

b <7 - v p 35 p 43 _ p 34 p 53 

3 4p 45 p 23 5p 54 p 23 

„ _ p 35 p 43 p 34 p 53 

3 " y4p 45 p 23 y5p 54 p 23 


(72) 


Note . In order to get these expressions for a 3 , b 3 , c 3 take into account 
the relation 


p ik " p ki 
and the identity 

p 23 p 45 +p 24 p 35 +p 25 p 43 = 0 

The three rows V^, W 2 , W 3 are dependent, since for 


A, = - 43 and A 7 = - 

1 23 2 23 


(73) 


(74) 


(75) 
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one has 


+ A 2 W 2 = W 3 (76 ^ 

Therefore when the five stations are coplanar the rank of the matrix 
of Table 8 is less than nine. That is, the plane determined by any three 

station-position points is a determinantai locus for the rest of the 
two station-position points. 

1.4.4 The Case of Six and More Ground Stations 

In the previous two cases, (i.e., in the cases of four and five 
ground stations) the discussion was more or less exhaustive because 
the number of basic patterns of observation was small (1 in the 
case of four stations, and 6 in the case of five stations). This will 
not be the case for six and more ground stations, for the number of pat- 
terns of observation is large (see Section 1.1.4, relation (7)). How- 
ever in order to find out the determinantai loci for any pattern of ob- 
servation in the case of six and more stations, one may follow the same 
routine work that was followed in the previous cases. That happens be- 
cause, as in the case of five stations, one can form the condensed jaco- 
bian matrix of any pattern of observation in the case of six or more 
stations by putting row-wise one after the other the condensed jacobian 
matrices of station-quadruples. For example the condensed jacobian ma- 
trix corresponding to the observational pattern when all six stations 
observe simultaneously a set of satellite-position points, is formed by 
putting row-wise one after the other (regardless order) the condensed 
jacobian matrices of any three station-quadruples which have three sta- 
tions in common. The set of the satellite-position points for each one 
of those three station-quadruples being the very same set observed by all 
six stations. Before proceeding further, a list of all possible basic 
observational patterns in the case of six stations, any 
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observational pattern consists o£, is given below for reference purposes, 
(see Section 1.1.4). 


List of the Basic Observational Patterns in the Case of Six Stations 


Quadruples 

Observed 

Fivefolds 

Observed 

Sixfolds 

* - -* ^ S*i 

Observed 

(15) 

Point-sets] 

(6) 

Point-sets 

CD 

Point-sets 

0 56 =(1, 2,3,4) 

<j S 6> 

F 6 =(l,2,3,4,5) 

0 6 } 

E= (1 ,2,3 ,4,5,6) 

0> 

0 46 =(1,2,3,5) 

046) 

F 5 =(l,2,3,4,6) 

05} 



0 45 =(1,2,3,6) 

O45 } 

F 4 =(l,2,3, 5,6) 

04} 



°36 = 0>2 » 4 > 5 ) 

0 3 6> 

F 3 e(1,2,4,5,6) 

03} 



°35=C1>2*4,6) 

0 35 } 

F 2 =(1,3,4,S,6) 




0 34 s(l, 2,5,6) 

O34 I 

F 1 =(2,3,4, 5,6) 

Oi> 



026=0.3,4,5) 

0*26 } 





025=(1> 3,4, 6) 

025> 





0 24 =(1,3,5,6) 

O2 4 J 





023 = O» 4,5,6) 

O23 } 





0 16 e(2,3,4,5) 

Ol6 J 





0 15 = (2, 3, 4, 6) 

Ol5 





0 14 =(2,3,5,6) 

0.14 J 





0 13 =(2,4,5,6) 

0'l3 \ 





0i2=(3,4, 5,6) 

0l2> 


1 


i 


For the solution of the problem of determining the coordinates of 
the six stations, the observational pattern must consist at least of two 
station-quadruples comprising all six stations, (say 0s6=(l,2,3,4) and 
0 3 4=(1, 2, 5,6) ) , and at most all possible ones, i.e., fifteen 
(C4 = 6!/4! 2! = 15) . 

Consider any ovservational pattern in the case of six stations, and 
let (Qj } be the set of all the satellite-position points associated with 
this pattern. If all six stations were observing simultaneously the set 
{Qj }, one would have the maximum number of observation equations possible 
with the set (Qj). That is, the set of the observation equations of any 
observational pattern in the case of six stations is a subset 
of the set of the observation equations obtained when 
all six stations observe simultaneously the same set of 
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satellite-position points. For a unique least-squares solution the rank 
of the condensed jacobian matrix under an admissible set of constraints 
must be twelve. Therefore, by considering the condensed jacobian matrix 
as if all stations observed simultaneously the totality of the satellite- 
position points, no matter what the actual observational pattern is, and 
checking its rank, if it is less than twelve, then this is the case for 
the rank of the condensed jacobian matrix of the actual observational 
pattern as well. However, if it is equal to twelve, the rank of the 
condensed jacobian matrix of the actual pattern may or may not be equal 
to twelve. Note that the rank of the condensed jacobian matrix for six 
stations can not be greater than twelve, for under an admissible set of 
constraints it has twelve columns. In other words, if a determinants 1 
locus of the observational pattern when all six stations observe simul- 
taneously a certain set of satellite-position points occurs, then that 
is also a determinantal locus of any other observational pattern with 
the same set of satellite-position points. As it was said earlier, the 
condensed jacobian matrix of the observational pattern when all six 
stations observe simultaneously a certain set {Qj} of satellite-position 
points is formed by putting row-wise one after the other (regardless of 
order) the condensed jacobian matrices of three station-quadruples having 
in common three stations and each observing the same set (Qj). 

From the above discussion, the reason why the observational 
pattern when three station-quadruples carry out all the observations 
is so important, becomes obvious. The determinantal loci of this 
observational pattern are searched in the remainder of the present section 

The symbols to be used are completely analogous to those previously 
used in the cases of four and five stations. Thus the following symbols 
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will be used: 

Aik& : For the determinant of fourth order, whose rows from the top 
to the bottom are the affine coordinates of the satellite-position point 
Qj = (Xj ,Yj , 1) and three of the six station-position points except 
P i ,P k* P l» * n orcler increasing index respectively. For example 



X, 

Y J 

2 - 

1 





.1 

3 





A 123= 

x 4 

Y 4 

z 4 

li 



(77) 


x 5 

ys 

z 5 

1 





X 6 

x 6 

Z 6 

1 




Also 

x i 

y± 

z i 

1 




mi,m 2 ,m3 
D i = 

% 

S 

\ 

1 



(78) 


Xm 2 

Y m 2 

z m 2 

1 





Xm 3 

Ym 3 

Zm 3 

1 




where 








i = 1,2, 3, 4 

,5,6 







and 








mj jn2> m 3 an y tri pl 

e of 

indices 

of 

the 

satellite-position points 

ob- 

served by statipn Pj. 







The condensed jacobian matrix 

of the case when three station-quad- 

ruples, comprising all 

six 

stations 

, carry 

out all the observations, 

is 

given by the following Tabl 

e 9, This case 

includes, as it was said 

above 

that when all six stations 

observe 

simultaneously all the satellite- 



position-points. Take the station-quadruples 

0 56 =( 1 , 2 , 3 , 4 ), 0 46 =( 1 , 2 , 3 , 4 , 5 ), 0 45 =( 1 , 2 , 3 , 6 ) 
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which observe the following point-sets respectively 
^ 56 } ’ ^46 ? * ^45 ^ 4 

If all six stations observe simultaneously all the satellite-position 
points, then 

^56 } E {j 46 } = ^ 45 . } ^ 

Adopt as in the previous cases the set of admissible constraints 
x l>yi» z l>y2> z 2> z 3 * constant 

For the determinantal loci the rank of the matrix of Table 9 is less than 
twelve. Thus one has to examine the zeroes of the determinants of all 
12x12 submatrices. It is repeated again here, that these determinants 
are invariants with respect to the affine transformations and rational 
integral functions of the coordinates, therefore they have to be express- 
ible as polynomials in the fundamental invariants (see Appendix F). In- 
deed one obtains by Laplace expansion, as in the previous cases, the fol- 
lowing expression for the determinant of a typical 12x12 submatrix under 
the above adopted set of constraints. 




|A(0 s 6.0 46 ,04 5 )|« 

P 1 CP 2 ^ 


Pl(P2(03(P4^) 


C-l)° l +o 2 +o 3 +0 4 


y 3 i 
( X W2 Ym 2 1 


^m 4 ,m 5 yn 6 ^m 7 .mg.mg ^WiO^ll *^12 


where 


mi ra 2 » 3 n 4 m s m 6 m 7 m 8 m 9 m 10 m u 
X A 2Ag A 3v^ A 3n6 i 4?6 (i 456 A 4S6 A 546 A 54 6 A 546 A 645 A 645 A 645 

«1 cU 5 6>U{j 4 6> u{ ^5 1 
»2» m 3 


e<j56> 
*7»"S^9 C U46* 


(S6), 


OuMvCM* «>• 


(46), 


C4S), 


if m 1 ,m 2 ,ni3 t{j S6 > 
if mi ,^2 c ( ^46 ) 
if mj,ni2*®3 


respectively, 

respectively. 

respectively. 


( 79 ) 


( 80 ) 
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Table 9 . The Condensed Jacobian Matrix for Six Stations when Three Station-Quadruples, 0^, an< ^ °45> 

Carry out all the Observations, 































1) Determinantal Loci for the Satellite-Position Points, 

Loci of second degree . Each 12 x 12 submatrix of the jacobian matrix 
must have at least three rows from each of the three station-quadruples 
and at most six rows. Working as in the case of five stations, one has 
in the general case when each of the three station-quadruples observes 
six and more satellite-position points the following. 

The second degree determinantal locus consists of three lock L 56 , 

L 46 and L 45 for the sets of the satellite-position points (j 5 6 >, ^46 
and{j 4 ^} respectively. The locus for the point-set {355 Hs a second 
degree surface passing through the points: origin, p i> p 2 »^ 3>^4 and an ) 

two points from {j 56 >- In other words L 56 is not a general second degree 
surface, but it is such that any seven of its points determine it com- 
pletely. The equation of L S6 is taken by equating to zero the determin- 
ant of any 12x12 submatrix which has three rows from the station-quadruple 
Og^=(l f 2, 3,4) , the determinant being considered as a functional with 
respect to the coordinates of one of the satellite-position points ob- 
served by 0 56 h( 1 , 2 , 3 , 4 ), (i.e., one from the set {j 56 >). The locus L 46 
for the point-set {j 46 > is a second degree surface passing through the 

points: origin, Pi,P 2 > p 3> p 5 and an X two P oints from { ^46 ThuS 146 

too, is not a general second degree surface, but it is such that any 
seven of its points determine it completely. The equation of L 4 6 is 
taken by equating to zero the determinant of any 12x12 submatrix which 
has three rows from the station-quadruple 0 46 =(1 , 2, 3 , 5 ) , the determinant 
considered as a functional with respect to the coordinates of one 
of the satellite-position points from the point-set {j 45 ) . The locus L 45 
for the point-set { j ' } is a second degree surface passing through the 
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points: origin Pi>P 2 > p 3’ p 6 an< ^ any two P oints f rom {345}- Therefore 

L45 is not a general second degree surface but it is such that any seven 
of its points determine it completely. The equation of L45 is taken by 
equating to zero the determinant of any 12x12 submatrix which has three 
rows from the station-quadruple 045^(1,2,3,6), the determinant being con- 
sidered as a functional with respect to the coordinates of one of the 
satellite-position points from the point-set (345). 

Note . It is stressed that the loci L 56 , L 46 , L 4 5 go together. That means 
that if 

the points {355} lie on 

the points ^45} lie on L46, and 

the points {345} lie on L45 

simultaneously , then the rank of the jacobian matrix is less than twelve. 
In other words the three loci L56, L 4 ^, and L45 constitute simultaneously 
a determinantal locus for the point-sets (j^), ^ 46 ^’ ^45^ respectively. 

In order to detect the existence of that determinantal locus one may 
determine each of the loci L55, L 4 6, L 4 5 as the surface which passes 
through nine of its points and verify that the same surface is determined 
by any seven of its points. Take for example the locus L55 for the point- 
set <W- This is a second degree surface which contains the points : 
origin, Pi>P2> p 3> p 4 and b in case of determinantal locus. The locus 
L56 in general has the form 

ailX^+a22Y^+ a 33Z^+2ai2^Y+2ai3XZ+2a23YZ+2ai 4 X+2a24Y+2a34Z+a 4 4 = 0 
and any nine of its point not in singular configuration are sufficient to 
determine the coefficients a^j . If the result of the estimation of ajj is 
such that two of them are identically zero, that means that any seven of 
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the points of L55 determine it completely. In case of determinantal locus 
this must be the case for all three Lc^, ^ 46 » an< ^ ^45 

Linear loci . For the linear loci one has to carry out the same routine 
work, as in the previous cases, that is to equate to zero each factor of 
the typical term in the expression ( 79 ). Thus one has 

(i) From 


x 2 1 

V 1 


°, He* j 45 }) 


( 81 ) 


it is concluded that the plane through the point P 2 = (x 2 , y2» z 2^ and perpen- 
dicular to the x-axis is a determinantal locus for all the satellite- 
position points, i.e., { ) U ( j UU45 }• 

(ii) From 


Y 3 

m2 Y m 2 


*3 

X 

X 


m3 Y m3 


= 0, (m2,m3e{j 56 }u(j46KHj 45 } 


( 82 ) 


it is concluded that each plane of the pencil of planes with axis the 
projection line of (^(xs.yj.zj) to the (x,y)-plane, is a determinantal 
locus for all the satellite-position points ({j 56 } U( j 46 ) U 0 45 }) • 

(iii) From 

^m 4 ,m5,m 6 _ ( m4f m 5 ,m 6 e { j $6 }) ( 83 ) 

it is concluded that each plane of sheaf of planes through the point 
(x 4 ,y 4 ,z 4 ) is a determinantal locus for the points {j^h 

(iv) From 

my _ e r . , (*0/1 \ 

D5 = 0 , (m7,m 8 ,m 9 e{j 46 }) ( 84 ) 

it is concluded that each plane of the sheaf of planes through the point 

p =(x 6 ,y 6 ,z & ), is a determinantal locus for the points {j 45 }. 
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( 85 ) 


(v) From 

j™10» n ll» m 12 = o, ( m 10 » m 11 J m 12 e{j 4 5 >) 
it is concluded that each plane of the sheaf of planes through the point 
(x^.y^, z&) , (is a determinantal locus for the points f 5 ) . 

(vi) From 

I m l E { j 5 6 1Ut:i 46 }U{j 45 } 

1(56) , if mie (j 56 } 

1 

Xu )=j(46), if mje {j 46 } (86) 

| (45) , if mje (j 45 } 

it is concluded that the planes (Pi,P3,P 4 ), (Pi>P3> p s) an< i ( p l> p 3> p 6) 
constitute simultaneously a determinantal locus for the point-sets (j^), 

{;i 46 } ’ ^45 } respectively. 

(vii) From either one of the relations 

ra2.m3 e (i 56 > U(j 46 } U(j 45 ) (87) 

1 (56), if m 2 ,m3 e {j } respectively 
(46), if m2, m3 e (j 46 ) respectively 
(45), if m2, m3 e { j 4 jj } respectively 
it is concluded that the planes (Pi,P 2 ,P 4 ), (Pi,P2> p s) and (Pj.P^P^) 

locus for the point-sets 

(m 4 ,m s ,m 6 e { j 56 }) (88) 

is a determinantal locus for 

(m 7 ,m 8 ,m 9 e {j }) (89) 


constitute simultaneously a determinantal 
(j 46 (j 45 > respectively. 

(viii) From either one of the relations 

A m 4 _ ri _ n A^6 

A 456 " °> A 456 - °> A 456 °» 

it is concluded that the plane (Pj,P2, p 3) 
the points 

(ix) From either one of the relations 
A 546 = °> a 546 = °> A 546 = °> 
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it is concluded that the plane (Pf»P2>^3) is a determinantal locus foT 
the points {j^}. 

(x) From either one of the relations 

A 645 = °» A 645 = °» ^45 = °» t m 10> m ll> m 12 e ^45* (90) 

it is concluded that the plane (P^,P2)P3) i- s a determinantal locus for 

the points { }. 

2) Determinantal Loci for the Station-Position Points 

As in the previous cases of four and five stations, one is interested 
in finding out determinantal loci for the station-position points indepen- 
dently from the satellite-position points. From the discussion of the 
determinantal loci for satellite-position points, it is concluded that 
there is no second degree determinantal locus for the station-position 
points independent from the satellite-position points. Therefore one 
proceeds to examine the linear ones, i.e., those consisting of straight 
lines or planes. 

The determinantal loci consisting of straight lines, are obtained in 
a way similar to that which was followed in the case of five stations. The gen- 
eral procedure which was described there, applies here although it is not 
needed, for the condensed jacobian matrix of the observational pattern 
under discussion is available. Thus, under the set of admissible con- 
straints x 1 ,yi,zi,y2,Z2* z 3 = constant, one has the following. 

(i) The straight line determined by Pj and P3 is a determinantal locus 

for P4,Ps and P&. 

(ii) The straight line determined by Pi and P2 is a determinantal locus 
for P4,P5 and P^. 

(iii) The straight line determined by any two of the three points Pi>P2’* > 3 
is a determinantal locus for the third one. 
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If at least one of the coordinates xj,yj,zj of P^ were unconstrained, 
then the straight line (P 21 P 3 ) would be also a determinantal locus for 
P 4 j P 5 and P 6 . 

Note: A general procedure of determining such a determinantal locus for 

any number of stations and any observational pattern, was given within the 
discussion of determinantal loci for station-position points in the case 
of five stations. 

For the determinantal loci consisting of planes one observes that, 
since in this observational pattern the three station-quadruples involved, 
(0^^=(1,2,3,4) , 0 4 £= (1 , 2 ,3, 5) , 0 4 j=(l , 2, 3, 6) have in common the three 
station-position points Pi,P2 and P 3 , if each of these quadruples is cop- 
lanar, then their planes coincide. Therefore one may start off searching 
for these loci by assuming all six station-position points are coplanar. 

Assuming that the six station-position points are coplanaT, consider 
the determinants 

- Ci.M' = 1,2, 3, 4,5, 6) 
expanded according to the first row, i.e., 

A iU = a iU x j + e ik£ Y j +Y ikt z j + (5 ik£ ( 91 ) 

where A ikP are functions of the station coordinates 

only, (compare with (33) and (68)). 

By putting 


222 1/2 
Pikj. = £ ( a ikJ£+ s ikt + TikS.) > 


(e = si B n of a{ R£ ) 


(92) 


one has 


pT^- = d, (for all i,k,£ and j) (92) 

where d is the distance of Qj i(Xj , Yj , Zj) from the plane of the station- 
position points. 
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Consider now the jacobian matrix of Table 9 along with the set of 
admissible constraints adopted above, (i.e., x i »yi > z 1 >Y2> z 2» z 3 =constant), 
and effect the following rank-equivalent operations. 

3 56 

(i) Divide each row j 5g of 0s6=(l, 2,3,4) by the corresponding A 456 

3 46 

(ii) Divide each row j 46 of 0 46 = (1, 2, 3, 5) by the corresponding A 546 

j 45 

(iii) Divide each row of 045 ^( 1 , 2 , 3 , 6 ) by the corresponding A ^ 4 5 

(iv) Subtract from X 2 -column, P256^ p 456 tiroes the x 4 ~column 

(v) Subtract from X 2 -column, p 246 //p 546 times the x 5 -column 

(vi) Subtract from X 2 -column, p 246/ p 645 times the X 6 ~column 

(vii) Add to xj-column, p 356/ p 456 times the X 4 ~column 

(viii) Add to y 3 ~column, p 356/ p 456 times the y 4 -column 

(ix) Add to x^-column, p 346/ p 546 times the x^-column 

(x) Add to y 3 -column, p 346/ p 546 times the y 5 -column 

(xi) Add to X 3 -column, 0345 / p 645 times the x 6 -column 

(xii) Add to y 3 -column, p 345 / p 645 times the y 6 -column 

(xiii) Subtract the first row in each station-quadruple from the rest 
rows in it. 

Then it is obtained the following rank equivalent matrix, (Table 10). 

The column-block of the first three columns of Table 10 is indepen- 


dent of the satellite-position points and has three non-zero rows, namely 



( 94 ) j 

w^[c* 2 -x 5 )^, CX5-X3 

( 94 ) 2 

M(x 2 -* 4 -|§> (/6-y 

(94)3 
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= ( 1 , 2 , 3 , 6 ) 0 ,, 5(1 , 2 , 3 , 5 ) 



Table 10 . The Transformed Condensed Jacobian Matrix for Six Coplanar Stations, when the 
Three Station-Quadruples 0 56 , 0^, 0^ Carry out all the Observations. 


















The determinant of the matrix of these rows, i.e., 

W x *3 X3 1 

det ' P 456 > 5 U 45 P2S6P346P345 (X x-X4) X S y 5 1 

W 3 [ x 6 ^6 1 

X 3 y 3 1 

x 4 y 4 1 ( 95 ) 

x 6 ^6 1 

x 3 y 3 1 

+ P 24S P 356 P 346 (x 2 -x 6 ) x 4 y 4 1 

x 5 X5 1 

is not identically zero. Therefore the plane determined by any three 
station-position points is not in general a determinantal locus for the 
rest three points. 

Since det [W^ ,W.?, W^] is invariant under affine transformations (see 
Appendix C) by equating it to zero one may interpret it geometrically. 

As a matter of fact by considering that determinant as a functional with 
respect to the coordinates of one of the station-position points, say 
x 6 ,y 6 ,z 6 , and equating it to zero, one has a second degree planar locus 
for P 6 5(x 6 ,y 6 ,z 6 ) . That locus passes through the points p i > p 2> p 3> p 4> p 5 
as it is verified immediately. If the plane of the station-position 



- p 246 p 356 p 345 (x 2 -x 5 ) 



points has some specific positions with respect to the coordinate system, 
the determinant (95) is zero no matter where the six points lie in it. 
Thus if the plane of the station-position points is perpendicular to the 
(x,y) -plane, whence the projections of the six station-position points 
on the x,y-plane are collinear, then the three determinants on the right 
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hand side of (95) vanish. Therefore the plane determined by any three 
station-position points is not a determinantal locus for the other points, 
unless in addition to coplanarity condition the determinant (95) vanishes. 
Remark. For any observational pattern in the case of six stations, con- 
sisting of at least three station -quadruples, the plane determined by any 
three station-position points, is not a determinantal locus for the other 
groups of three. 

Suppose now that one has the case of n stations. From the procedure 
of getting the condensed jacobian matrix, it is concluded that the con- 
densed jacobian matrix of the observational pattern when all n stations 
observe simultaneously a certain set (Qj } of satellite-position points, 
is formed by putting row wise one after the other (regardless order) the 
condensed jacobian matrices of (n-3) station-quadruples which have in 
common three stations, and each one of which observes the very same set 
(Qj } of satellite-position points. Therefore the condensed jacobian 
matrix of any observational pattern (for any n>4) is formed by putting 
row wise one after the other, the condensed jacobian matrices of station- 
quadruples, and consequently for the determinantal loci one may repeat 
the same routine work that was followed in the previously treated cases. 

From the procedure of finding out the determinantal loci of an obser- 
vational pattern which was established in the previous discussions via 
forming the condensed jacobian matrix, one may describe general procedures 
of finding the determinantal loci of an observational pattern, without 
having to form the condensed jacobian matrix. These general proce- 
dures are described in the following section. 
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1.5 Synopses 


Let 

{P i }-(P 1 ,P 2 ,...,F > n ), ie{i}s(l, 2, . . . ,n) 
be a set of ground stations, observing in some observational pattern a 
certain set 

(Qj } = CQl » ■ • • j Q m ) > jet j ) = (1> 2 . ,m) 

of satellite-position points. 

I. Station Quadruples Composing the Observational Pattern 

Each satellite-position point Qj is observed by A stations, where 
4SA<n. 

The observation of Qj by X stations is equivalent to the observation of 
Qj by (A-3) station-quadruples formed from the X stations and having in 
common three stations. Note that there are 

C X = X! = X(X-l) (A-2) 

3 (A-3) ! 3! 6 

selections of (A-3) station-quadruples formed from X stations and having 
in common three stations. Therefore the set of station-quadruples com- 
posing an observational pattern is not unique as long as there are satel 
lite-position points observed by more than four stations. 

For the total number of the station-quadruples composing an observa 
tional pattern one has to group the totallity of the satellite-position 
points in groups, each one of which is observed by the same stations. 
Thus, suppose that 

(Qj h> {Qj h* ■ • { Qj 

be sets of satellite-position points observed by 

^1 > ^2> ■ • • i 
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stations respectively. It is understood that 

{Qj Ji { Qj } 2 *'* {Qj }k=tQj 1= (Ql»Q2. • • • iQm) • 

Then the totality of station-quadruples for that observational pattern is 

k 

K = CX r 3)+(X 2 -3) + --* + CX k -3) = jijtti-3) 
or k 

K = Z X\ -3k, 
i=l 

For example, if there are only four stations, then all the satellite- 
position points are observed by the same four stations and in this case 
k=l, that is 

{Q j }lE{Q j }=(Q 1 ,Q 2 ,...,Q m ). 

The formula for K above gives 
K = 4-3 = 1. 

II. The Condensed Jacobian Matrix of a Station-Quadruple 

This is given in Table 1 (see Section 1.4.2). The same table is re- 
called here with a change in the indices of the determinants . Thus 
instead of the subindex i which indicates the station whose coordinates 
do not appear in the determinant, one uses the indices of the stations 
whose coordinates appear in the determinant. Without loss of generality 
take the s tat ion -quadruple Pi,F > 2*f , 3>f , 4’ Then the following identities 
exist (see relations [27]). 


4=434,4=^34, A I =A 1 24, = A 1 2 3 , 


and Table 1 is written 
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Notes. (1) The condensed jacobian matrix of a station-quadruple has four 
column-blocks, one for each station, with three columns each block corres- 
ponding to the x,y, 2-coordinates of the respective station. (2) All 
three-column blocks have the same form. Thus in general the three- 
column block corresponding to station P^ in the condensed jacobian matrix 
of a station-quadruple including Pj, say P] c ,P£,Pi,P m ,(k<«.<i<m) is given 
by 

e (Xj -Xi) e (Yj -y±) £ -z i) > 

where e = ±1, according to the following stipulation. 

! fl, if i is first or third in order among i,k,£,m 

. 

-1, if i is second or fourth in order among i,k,fc,m 
(3) The condensed jacobian matrix of a station-quadruple has as many 
rows as satellite-position points observed by that quadruple. 

Ill . The Condensed Jacobian Matrix of any Observational Pattern 

The condensed jacobian matrix of any observational pattern is formed 
by putting row wise one after the other (regardless order) the condensed 
jacobian matrices of the station-quadruples composing that observational 
pattern. 

A symbolic scheme of the condensed jacobian matrix is given below, 
which is useful for the description of the general procedures of deter- 
mining the determinantal loci without having to write down the explicit 
form of the condensed jacobian matrix. That is established through the 
following three examples. 

Example 1 . Let an observational pattern in the case of five stations 
(Pi , P2, P3, P4) be composed by the station-quadruples 
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f 


°1234 S ( P 1> P 2’ P 3’ P 4)» and 0 12 35=(P 1 ,P2. p 3’ P 5^ 

Then the condensed jacobian matrix may be written in the following sym- 
bolic way. 


- — 
°1234 


p l 

P 2 P 3 P 4 o' 

°1235 


p l 

p 2 p 3 0 P 5 


Example 2 . Let an observational pattern in the case of six stations 
( p l » p 2» p 3> P 4 ,P 5* P 6^ b e composed by the station-quadruples 

°1234=( p l> p 2> p 3> p 4)» °1256 E ( P 1> P 2* P S* P 6) » °3456=( P 3> P 4> P 5* P 6) • 

Then the condensed jacobian matrix may be written in the following sym- 
bolic way. 


°1234 


’ p l 

p 2 

p 3 

p 4 

0 

0 " 

°1256 

= 

p l 

p 2 

0 

0 

p 5 

p 6 

°3456 


0 

0 

p 3 

p 4 

p 5 

p 6 


Example 3. Let an observational pattern in the case of six stations be 
composed by the station-quadruples 

°1234=( p l> p 2> p 3* p 4) » 0 1235=( p l’ p 2> p 3> p 5) 

°1256 i ( P l> P 2> P 5^ P 6)» °2345-( P 2* P 3’ P 4> P s)* 

Then the symbolic scheme of the condensed jacobian matrix is 


°1234 


p l 

p 2 

p 3 

p 4 

0 

0 

°1235 


p l 

P 2 

P 3 

0 

P 5 

0 

°12S6 


p l 

P 2 

0 

0 

P 5 

P 6 

°2345 


0 

_ 

p 2 

p 3 

P 4 

p 5 

0 


IV. Definition of the Determinantal Loci 

Consider an observational pattern involving n stations. Its con 
densed jacobian matrix has n three-column blocks, i.e., 3n columns. 
The rank of this matrix is at most 
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y = 3n-6 = 3 (n-2) . 

Therefore for a unique least squares solution of the problem, one has to 
introduce six relations independent from each other and the rows of the 
condensed jacobian matrix- In the previous discussion such relations 
were introduced as follows: 

x l»yi» z 1^2^ z 2» z 3 = constant - 

A determinants 1 locus is that relation between the coordinates of 
the satellite and/or station position points for which the rank of the 
condensed jacobian matrix is less than y = 3n-6 = 3 (n-2), or that is the 
same, that which annihilates the determinants of all (3n-6)x(3n-6) sub- 
matrices. 

Remark: For the determinantal loci one has to search for the zeroes of 

the determinants of all (3n-6)x(3n-6) submatrices. Without loss of gen- 
erality adopt the admissible set of constraints. 
x l’^l ,z l , ^2 jZ 2» z 3 = constant. 

Then in order the determinant of a (3n-6)x(3n-6) submatrix not to be 
identically zero, it must have from each station-quadruple a minimum num- 
ber of rows. That minimum number of rows, depends upon how many stations 
out of the four in that quadruple do not participate in any other station 
quadruple of the same observational pattern , and how many coordinates of 
those stations are unconstrained. It varies from six , when two uncon- 
strained stations or three stations whose three appropriate coordinates 
are constrained, or four stations whose six appropriate coordinates are 
constrained, do not participate in any other station-quadruple of the 
same observational pattern, to one when all four stations of the station- 
quadruple in question participate in at least one more station quadruple 
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of the same observational pattern. This remark is of importance for the 
second degree determinantal loci for the satellite-position points, (see 
the following synopsis V). 

V . Description of the Determinantal Loci 

The determinantal loci are classified as follows. 

1. Determinantal loci for the satellite-position points 

la. Determinantal loci of second degree.. 

lb. Linear determinantal loci. 

2. Determinantal loci for the station-position points. 

2a. Determinantal loci consisting of planes. 

2b. Determinantal loci consisting of straight lines. 

Itemark : The determinantal loci for the station-position points are inde- 

pendent from the satellite-position points, and so one does not have sec- 
ond degree determinantal loci for the station-position points.* 

In the following, each of the above classes of determinantal loci 
is described in a general way for any number of stations and any observa- 
tional pattern, without utilizing the explicit expression of the condensed 
jacobian matrix, but only the symbolic one, (see synopsis III above). 

1 . Determinantal loci for the satellite-position points , 
la. Determinantal loci of second degree . 

Let an observational pattern involving n stations 
P 1 ,P 2 ,...,P n , (ns4) 

a set of m satellite-position points. 

{Qj} = (m>[3(n-2)/(n-3)j),. 

and being composed by a certain set of station-quadruples, 

*The second degree determinantal loci involve both satellite and station- 
position points. 
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[n/4]*+a <p<n! / (n-4) ! 4! 

1 0, if n/4 is an integer 
1, if n/4 is not an integer 

Let also 

<Vl’ 1Q j } 2 'Vp 

be the sets of the satellite-position points observed by the station- 
quadruples respectively. It is understood that these sets con- 

stitute a decomposition of the set {Qj}, i.e., 

{Qj }] U (Qj }2 U * ’ * U(Qj}p= {Qj } = (QiiQ2» • • ■ »Qj > • • • » ■ 

Then one may determine a set of p second degree surfaces (some of which 

or even all of them may coincide), 

Si,s 2 ,--.,s p 

one for each station-quadruple, such that they constitute simultaneously 
a determinantal locus for the sets of the satellite-position points 
{Qj {Qj • • * » * p 
respectively. That is, if 

{Qj U £S 1 * {Qj^2^"^2’ {Qj } p eS p 

then the determinants of all (3n-6)x(3n-6) submatrices of the condensed 
jacobian matrix vanish, and consequently the rank of that matrix is less 
than 3n-6. 

One now proceeds to describe the procedure of determining the second 
degree surface 

S v , (v=l,2, . . . ,p) 

corresponding to the station-quadruple 0 V . For some fixed value of v let 
0y=0q = (Pq , Pj. , Pg , P-^ ) , q,r,s,te (1,2, . . . ,n) 

*The square brackets indicate the integer part of the bracketed number. 
tThe symbol e stands for expression "belongs to". 
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The determination of the surface S v , depends upon the minimum number of 
rows allowed by the observational pattern to be taken from the quadruple 
° v E 0qrst whlle forming the (3n-6)x (3n-6) submatrices, in order that their 
determinants not be identically zero. That minimum number of rows 
depends upon how many stations out of the four in that station- quadruple, 
do not participate in any other station-quadruple of the same observa- 
tional pattern, and how many coordinates, if any, of those stations are 
unconstrained. It varies from six to one as follows: 

It is six, when (a) two stations of the quadruple without 
constrained coordinates, or (b) when three stations of the quadruple 
with three appropriate* coordinates constrained, or (c) when all four 
stations^ of the quadruple with six appropriate coordinates constrained, 
do not participate .in any other station-quadruple of the same obser- 
vational pattern. 

It is five, when two stations of the quadruple, one with uncon- 
strained coordinates, and the other with one out of the three coordinates 
constrained, do not participate in any other st at ion -quadruple of the 
same observational pattern. It is four , when two stations of the quad- 
ruple, one with unconstrained coordinates, and the other with two out of 
the three coordinates constrained, do not participate in any other 
station-quadruple of the same observational pattern. 

It is three , when just one station of the quadruple, with uncon- 
strained coordinates, does not participate in any other station-quadruple 
of the same observational pattern. 

*By appropriate one means that these constraints are members of an ad- 
missible set of constraints. 

tin that case the quadruple is independent from the others if there 
exists any. 
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It is two , when just one station of the quadruple with one out of its 
three coordinates constrained, does not participate in any other station- 
quadruple of the same observational pattern. 

It is one , when one station of the quadruple with two out of its 
three coordinates constrained, does not participate in any other 
station-quadruple of the same observational pattern. 

In case all four stations participate in other station-quadruples 
of the same observational pattern, one may form (3n-6)x(3n-6) submatrices 
which do not include any row from that station-quadruple. That is why 
the minimum number of rows from that quadruple is one, even in this 
case . 

After the description of determining the minimum number of rows al- 
lowed by the observational pattern to be taken from a station-quadruple 
0 V , while forming the (3n-6)x(3n-6) submatrices, in order their determin- 
ants not to be identically zero, one proceeds in the description of de- 
termining the second degree surface S^, as the locus of the satellite 
position points 

{Qj 

observed by the station-quadruple 
^v z ^qrst “ 

The surface S v passes through the points: origin, Pq,P r ,Ps>Pt anc * the 

satellite-position points observed by the station-quadruple 0v=0qrst* Let 
t be the minimum number of rows from the station-quadruple 0 V . Then the 
surface S v is such that it is determined by the points: origin, 
Pq,P r ,P s ,P 1: and any (t- 1) satellite-position points. Since a general 
second degree surface entails nine of its points to be determined, it is 
concluded the following. 
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(i) Fort =6 or 5, S is a general second degree surface, 

(ii) For.T =4, 3, 2,1, S is not_ a general second degree surface, but 
it is such that any eight, seven, six, five of its points respectively 
determine it completely. 

Remark: The determination of a second degree surface by nine of its 

points, is a very well known problem in analytic geometry. One has to deter 
mine the coefficients of the analytic expression of a general second sur- 
face, i.e., 

a 11 X 2 +a22Y 2+ a33Z 2 +2a 12 XY+2a 1 3XZ+2a23YZ+2a 1 4X+2a24Y+2a34Z+a 44 = 0. 

If the second degree surface is not a general one, then some of the coef- 
ficients in the above expression are identically equal to zero. 

Thus the description of the procedure of determining the second de- 
gree determinantal locus is complete. Note that for a given observation- 
al pattern there exists just one second degree determinantal locus, con- 
sisting of the second degree surfaces , S 2 ,. . . ,S V . 

lb.. Linear det erminantal loci 

(a) Each plane of the sheaf of planes through a station Pi, with 
none of the coordinates is constrained, is a determinantal locus for the 
satellite-position points observed by all the station-quadruples where 
the station Pi participates in. If one out of the three coordinates of 
Pi is constrained, say Zi, then instead of the planes of the sheaf through 
p. one has the planes of the pencil whose axis passes through P^^ and is 
perpendicular to the (x,y) -plane. 

If two out of the three coordinates of Pi are constrained, for example 
y it z^, then instead of the sheaf of planes through P^, one has the plane 
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through and perpendicular to the x-axis. 

(3) Another category of linear determinantal loci for satellite- 
position points consists of sets of planes determined by the station- 
position points. Each determinantal locus of this category is associated 
with a station. In other words there are as many determinantal loci of 
this category as stations. Consider any station of the network, say P^. 
The determinantal locus associated with this station is determined as 
follows. 

Let 

°1»°2> • • • >°k 

be the station-quadruples where P^ participates in, and 

f Qj 5 1 > (Qj )2 > • • ♦ » ^ Qj ^k 

be the sets of the satellite-position points observed by those station- 
quadruples respectively. Excluding station-position point P^, the re- 
maining three station-position points in each one of the station-quadru- 
ples 0j,0 2 , ... , 0^ determine a plane. Let these planes be denoted by itj, 
ir 2 , . . . ,ir k . Then the planes irj ,tt 2 , • • constitute simultaneously a deter- 
minantal locus for the satellite-position point-sets (Qj )]_, (Q^ ) 2 , - ■ • » tQj 
respectively. That is, if 
{Qjll e tt x 

{Qj } 2 e it 2 

{ Vk £ *k 

there is no unique least squares solution. 

Remark : In case of a determinantal locus of this category, say that 

associated with station P^, the three columns of the condensed jacobian 
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matrix corresponding to the coordinates xp,yi,2p of become zero ones. 

These are the determinantal loci for the satellite-position points. 
One might add to the above quoted determinantal loci, the totality of 
the straight lines in space. However, these are included in the above 
quoted case (a), for then there exists always a plane through the line 
in space and a station-position point. 

2. Determinantal loci for the station-position points . 

2a. Determinantal loci consisting of planes . 

It has been demonstrated in previous discussions that in the case 
of four stations, the plane determined by any three out of the four 
station-position points is a determinantal locus for the fourth one. 

That is, if the four station-position points are coplanar there is no 
unique least squares solution. 

Also in the case of five stations, the plane determined by any 
three out of the five station-position points is a determinantal locus 
for the other two. That is, if the five station-position points are cop- 
lanar there is no unique least squares solution. 

In the case of six stations it has been demonstrated that, if the 
observational pattern is composed of at least three station-quadruples 
there exists in general no determinantal locus consisting of planes. In 
the specific observational pattern composed by three station-quadruples 
having in common three stations, the plane determined by any three out of 
the six station-position points is a determinantal locus for the other 
three station-position points, provided that plane is perpendicular to 
the (x,y) -plane, under the admissible set of constraints: 

X 1 »Y\ » Z 1 z 2' z 3 = cons tant. A remarkable note may be drawn from this 
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specific case. It is observed that the selection of the coordinate sys- 
tem may create a determinantal locus, which otherwise would not exist. 
Thus if in the case of six coplanar station-position points, one selects 
a coordinate system with the (x,y)-plane perpendicular to the plane of 
the station-position points, then the observational pattern composed by 
three station-quadruples with three stations in common, and under the set 
of constraints: Xj .yj , Zj ,y 2 , z 2 > z 3 = constant, does not yield a unique 

least squares solution. In general this observational pattern [involving 
six coplanar station-position points and being composed by three station- 
quadruples with three stations in common), does not yield a unique least 
squares solution, if in addition to the coplanarity condition, the gener- 
al second degree planar curve, determined by any five out of the six 
station-position points, passes through the sixth one, (see Section 
1.4.4). 

It should be noted that even in the case of six stations, one may 
encounter observational patterns where do exist determinantal loci for 
the station-position points consisting of planes. These observational 
patterns are composed by two station-quadruples. Consider for example 
the observational pattern whose the condensed jacobian matrix has the 
following symbolic scheme. 


r — 


- 

“ 

°1234 

— 

P 1 

P 2 P 5 p 4 0 0 

°1256 


Pi 

P 2 0 0 P 5 P 6 


For this observational pattern the plane determined by any three out of 
the four station-position points of either station-quadruple ( 0 12 3 4 and/or 
^1256^ * s a determinantal locus for the fourth one. 
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In a general observational pattern, the plane determined by any 
three out of the four station-position points of a station-quadruple, is 
a determinantal locus for the fourth one, if two stations of that quad- 
ruple with unconstrained coordinates, or three stations with three appro- 
priate coordinates constrained, or four stations with six appropriate 
coordinates constrained, do not participate in any other station quad- 
ruple of the same observational pattern. 

In the observational pattern involving n stations (n>6) and being 
composed by v station-quadruples (v>3) which have in common three sta- 
tions, provided the station-position points are coplanar, the second de- 
gree planar curve determined by any five station-position points, is a 
determinantal locus for the other (n-5) ones, 

2b. Determinantal loci consisting of straight lines 
Each of these determinantal loci is associated with a station. In 
other words there are as many determinantal loci of this type as stations. 
The general procedure of finding out such a determinantal locus is ana- 
logous to that which was described above for the satellite-position 
points. 

Consider any station of the network, say Pj. The determinantal 
locus for station-position points associated with this station is deter- 
mined as follows . 

Let 

0l> * • • * 

be the station-quadruples in which P i participates. By excluding P^, 
there remains from each one of the station-quadruples Oj , 02 , ♦ . * , 0^ a 
point-triple. Let these triples be denoted by 
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Ill 222 
{ p l> p 2> p 3*> ( p l> P 2> p 3^» 


r ^ k 

.... { p i, p 2 » 



Let also 

(ei), O2), •••> C £ k) 

be the straight lines determined by any two of the three points of the 
above point-triples respectively. Without loss of generality one may 
take 

(c j) = (Pj > » C e 2^ = ^ P 1 ’ P 2^ ’ ' " ’ C e k ) = (P l’ P 2^ 

Then the straight lines (e^), (£2), (^k) constitute simultaneously 

1 k 

a determinantal locus for the station-position points P3, . . . , P3 respec- 
tively. That is , if each one of the point-triples 

< P 1» P L P 3>> { p l. p 2> p 3>> •••» { p i, p 2. p 3> 
is collinear, the x^,y^ ,z^ -column block of the condensed jacobian matrix 

becomes a zero one, and there is no unique least squares solution. 



2. DETERMINANTAL LOCI FOR RANGE-DIFFERENCE NETWORKS 
2.1 Introductory Remarks 


The present discussion refers to the "continuously" measured range- 
differences. This observational mode is that of geoceiver - a new com- 
pact and portable device for satellite tracking. 

Let Qj, (j=0,l,2,...,m) be a set of satellite-position points at the 
instances t . , (j=0, 1 ,2 , . . . ,m) , and P-, (i=l, 2, 3,4, . . . ,n) be a set of 

J X 

ground stations observing by means of geoceiver. 


Q 0 Q, Q 2 Q 3 Qj-i Qj 



Let also sj be the range from to Qj , and the range difference be- 

i k 

tween the ranges and S^, i.e., 

s{ ,k = Si - s{ , (k>j) ( 96 ) 
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When station observes by means of geoceiver the satellite-position 
points Q 0 , Q r Q 7 , • Qj , • ^ (see Figure 3) then the observed 
quantities are the range differences: 


4- 1 


s ; 


1,2 c2 


s i 



gm,m-l 

i 


It is noted that these symbols stand for the true values of the denoted 
quantities. Because of the "continuity” one may find the range differ- 
ences of one range, say (PiQo)-Si from each one of the others. Thus one 
has 


s) 

1 


s9 + s?* 1 

i i 



+ Sj- 1 * J 


s9 + z 

1 k=l 1 


or 



q k-l,k 


(j — 1,2, .. ., m)% 


(j = l,2,3,. . . ,m) 

(97) 

(98) 


The right hand part of (98) is known from the observations. One may sub- 
stitute the set 


r <- ? > J cM c?, 3 cj l,j c 1 ? 1 

J j = l,2,...m - lc> i * 

(99) 

by the equivalent one 

tSi ,:l } jsl> 2,...,m = (S?- 1 , S?’ 2 , S-’ 3 , ..., S?’\ .... S?’ m ). (100) 


Indeed the set (99) is equivalent to (100) because if the members of the 
set (99), (being considered as functional of the coordinates) are alge- 
braiclly independent, this is also the case for the members of the set 
(100), for each member of (99) is a linear combination of members of (100). 
As a matter of fact 

= s9»j - sPJ- 1 (ioi) 
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In the following the set {S^ 3 ^ ^} = i t 2 , . ,m cons idered . 
Qj ~(Xj,Yj,Zj), j=0,l,2,...,m 
P i = fx i >yi> z O> 

Then 


s i >:} 58 s l ~ s i 

1/2 _ 0 7 1/2 
= [ (Xj -x ± ) 2 + (Yj -y t ) 2 + (Z^ -z ± ) 2 ] - [(X 0 -x i ) 2 + (Y 0 -y i ) 2 +CZ 0 -z i ) ] 

( 102 ) 


Since at each event three unknown coordinates are introduced (i.e., the 
coordinates of the satellite-position point being observed), at least four 
stations must observe simultaneously*. 

2.2 The Case of Four Ground Stations 
Let Pi=(xi,yi,Zi) , (i=l,2,3,4) be four stations observing simultane- 
ously by means of geoceiver the satellite-position points Qj , (j-1,2, . . ,m) . 
Then one has a set of 4m simultaneous equations of the form (102) , with 
3(m+4) unknowns, six of which are to be constrained in order to establish 
the coordinate system, so that there actually are 3(m+2) unknowns. 

For a unique least squares solution the rank of the jacobian matrix 
J M of the 4m functions (102) must be equal to 3(m+2), i.e., 

rank = u = 3(m+2). (103) 

As in the case of range networks, one may condensef the jacobian 
matrix. The condensed jacobian matrix will be obtained here via elimina- 
tion (see Appendix E) . Take the differential form of (102). 


*This must be the case for geometric satellite geodesy. 

+The term condensed jacobian matrix is used conventionally throughout 
this work, and it is not a general mathematical concept. It is used here 
on the grounds that the determinants of the submatrices for the determin- 
antal loci (not of any square submatrix) are equal to the corresponding 
ones of the condensed jacobian matrix. 
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dS?’ gj[ (Xj -x^dXj + (Yj -yjJdYj + (Zj - z ± ) dZ^ ] 
i 

- [ (h: x l _ x 9.- Xi )dx i+ ( Y r yi - Y °" yi )d yi 4 (hzli .. z °~ zi )d Zi ] 


S 1 


so 


S 1 


so 


S 1 


SO 


or 


;0, j_ 


dS i * a ij dX j +b ij dY j +c ij dZ j-( a ij- a i0^ dx I^ b ij- b i0^ d yi^ c ij- c i0^ dz i ( 104) 
Note that dXQ,dYg,dZQ= 0, for Xq,Yq,Zq may be chosen arbitrarily. How- 
ever the choice of XqjYq.Zo is going to be counted for the constraints in 
order to establish the coordinate system. 

Relations (104) render the mapping of the differential space spanned 

by dXj >dYj , dZj , dx^,dy^,dz^, (i=l,2,3,4; j-i, . . . ,m) into the differential 

0, j 

space spanned by dS^ , the matrix of the transformation being the jaco- 
bian matrix Jj^. Being interested in the rank of the matrix of a linear 
transformation, it is sufficient to consider the kernel of that transfor- 
mation. 

Thus one takes 

a ij dX j +b ij dY j +c ii dz ) ~ ( a ij -a i0^ dx i‘ ( b i j -b io) d Yi _ C c ij ~ c io) dz i = 0 ( 10S ) 

where 


a ij> b ij> c ij _ the direction cosines of Pj^Qj 
a i0» b i0 jC i0 _ tbe direction cosines of PqQq. 

By grouping the 4m equations (105) in m groups, one group per value of 
the index j, and eliminating from each group the differentials dXj,dYj, 
dZj, one comes up with the resultant system of m simultaneous equations, 
the matrix of the coefficients of which is what has been called condensed 
jacobian matrix. That is given in Table 12, where 


^1* ( i= l> 2 > 3 > 4 ; j = l,2,..,,m) 
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dx l 

d/i 

dz l ' 

dx 2 

*2 

dl 2 

dx 3 

mm 

dz 3 

dx 4 

dy 4 

dz 4 

( a ir a io)4 

1 ( b lj- b lo) A l 

(«ir c io)4 

-(a2j-a2o)4 

- 

-(b2j-b20^J 


Ca 3 j-a 3 o)4 


(Cjj * c 30)4 

-(a4j-«4o) 4 

-Cb4j-b40)4 

-( C 4j-C40) A J 

1 


Table 12. The Condensed Jacobian Matrix for Four Stations "Continuously" Observing Range Differences. 















stands for the same determinant as in the case of range networks. To 
recall it, A^ is given 


x i 

Y ; 

Z i 

1 

X 1 

Yl 

Z 1 

1 

x 2 

x 2 

z 2 

1 

x 4 

x 4 

z 4 

1 


One notices immediately that the entries of the condensed jacobian 
matrix are not rational functions of the coordinates of the satellite and 
station position points. Therefore one can not handle this matrix as 
that of range networks. To find out the determinantal loci in this case 
one may proceed as follows. First adopt a set of admissible constraints. 
Since the constraints Xq.YqjZq = constant have already been adopted, 
three more are needed. By constraining the coordinates yj,z^,Z 2 one has 
the following set of admissible constraints. 

x 0> y 0* Z 0^1’ z 1> z 2 = constant. (106) 

By leaving out from Table 12 the columns corresponding to dyj,dzj,dz2 and 
taking the determinant of a typical 9x9 submatrix, it is obtained by 
Laplace expansion the following expression (compare with (29). 


M - ]£ 

PlV>2(P3)) 




a 20 b 20 1 

a 10 * 



a ij . i 


a2 h b2 i: 1 

J 1 


82 J3 b2 j3 ' 


a 30 

b 30 

c 30 

1 

* 3 U 

b3 M 


1 

* S J 5 

bj JS 

C, i5 

1 

' S i6 

b3 J6 

CJ j6 

1 


a 40 

b 40 

c 40 

l 

a “i7 

b4 J7 

C4 37 

1 

34 j 8 

"to 

C4 38 

1 

84 39 

b< J9 

C< i9 

1 


X 


4 1 *2 *3 4 3 4 3 *4 4 4 % 


(107) 
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where 


j 1 J 2 ’ ,: *3 

is any selection of nine satellite-position points out of the m. 

The determinants 1 loci for the satellite and station position points 
are obtained from the expression (107) by equating to zero each factor of 
the typical term. Thus one has the following: 

(i) From 

11 


a 10 

a iii 


= o 


(108) 


it is concluded /that if all the satellite-position points, (Qo»Qi > • ♦ • » 
lie on the circular cone with vertex at Pi=( x i>yij z l) an ^ axis parallel to 
the x-axis, (when ajo^ljj* (j 1 = 1»2 . ,m)) , then the determinants of all 

9x9 submatrices vanish. Therefore the circular cone with vertex at 
Pl=(xi,yi, z i), axis parallel to the x-axis and passing through the point 
Qq= (Xq, Yq, Zq) is a determinantal locus for the satellite-position points 


CQl.02. ' * 4 4 

(ii) From 

| a 20 b 20 


a 2j ; 

a 2j, 


b 2j. 


= 0 


(109) 


’3 ” J 3 

it is concluded that if all the satellite-position points are distributed 
on any two semi-straight lines through P 2 =( x 2,y2 » z 2) » (which may coincide), 
then the determinants of all 9x9 submatrices vanish. Therefore any two 
semi-straight lines through P 2 E ( x 2 *Y2> z 2^ » constitute a determinantal 
locus for the satellite-position points Qq.Qj .... ,Q m . Since the two semi- 
straight lines may coincide, any single semi-straight line through P 2 » is 
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also a determinantal locus for the satellite-position points Qo,Qi , • - • ,Q m 

(iii) From 

a 30 b 30 c 30 

a 3j 4 b 3j 4 c 3.i 4 

a 3j s b 3j s =3j s 

a % b 3 j 6 cjj 6 

it is concluded that if all the satellite-position points are distributed 
on three semi-straight lines through P 2 E ( x 2 >/ 2 > z 2 ) » ( an >' two °f which or 
all three may coincide) then the determinants of all 9x9 submatrices van- 
ish. Therefore any three semi-straight lines through P 3 = (* 3 »y 3 > 23 ) con- 
stitute a determinantal locus for the satellite-position points 
Q 0 ,Qi,Q 2 ,...,Qm* Since any two of the three semi-straight lines or all 
three may coincide, any two or any single semi-straight line through P 3 , 
constitutes also a determinantal locus for the satellite-position points 

Qo»Qi>Q 2 » ♦ • • »Qm‘ 

(iv) From 

a 40 b 40 c 40 1 

a 4j 7 Vj 7 C 4j 7 1 

a 4j g b 4j g C 4j g 1 

a 4j 9 b 4j 9 C 4j g 1 

it is concluded as above that, any three semi-straight lines through 
P 4 = (X 4 ,y 4 , Z 4 ) , (any two of which or all three may coincide), constitute a 
determinantal locus for the satellite-position points Qo»Ql» • • • >Qm* Since 
any two of the three semi -straight lines or all three may coincide, any 
two or any single semi-straight line through P 4 , constitutes also a deter- 
minantal locus for the satellite-position points Qq,Qj, . . . ,Q m . 
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( 112 ) 


(v) From 

Aj 1 = 0 

it is concluded that if all the satellite position points lie in the plane 


(P 7 ,P 3 ,P 4 ), the condensed jacobian matrix is of rank less than nine and 
therefore the plane (P 2 > p 3 > p 4 ) is -a -determinantal locus for the set of all 
the satellite-position points. ■ 


(vi) From either one of the relations 

j 2 ^3 

^2 = 0 , A2 =0 


(113) 


it is concluded that the plane ( p i, p 3 , p 4 ) is a determinantal locus for the 
set of all the satellite-position points. 

(vii) From either one of the relations 

4 4 = 0 . 4 s - 0 , 4 6 -.0 ' ( 1 M ) 


it is concluded that the plane ( p i, p 2> p 3) is a determinantal locus for the 
set of all the satellite-position points. 

(viii) From either one of the relations 

4 7 = o, 4 8 - 0, 4 9 » o (115) 

it is concluded that the plane ( p i> p 2> p 3) is a determinantal locus for the 
set of all the satellite-position points. 

For the station-position points one has the following. , 

By considering the determinants a{, i=l, 2,3,4 one concludes that, if 
any three station-position points are collinear, the last three rows in 
one of these determinants are linearly dependent, and the determinant of 
any 9x9 submatrix of the condensed jacobian matrix vanishes. In other 
words the straight line determined by any two station-position points is 
a determinantal locus for each one of the other two separately. These are 
the only determinantal loci for the station-position points independently 
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from the satellite-position points. The plane determined by any three 
station-position points, is not a determinantal locus for the fourth one. 

2.3 The Case of More Than Four Ground Stations 
When more than four stations participate in the observations, one 
may consider the condensed jacobian matrix being formed from the con- 
densed jacobian matrices of station-quadruples placed row wise one after 
the other, in an exactly analogous way as in the case of range networks. 
Then the determinantal loci for the satellite-position points result by 
repeating the above quoted for each station quadruple envolved, and the 
determinantal lock for the station-position points too, result in the 
same way as in the case of four stations. Thus, the straight line deter- 
mined by any two of the station -position points, is a determinantal locus 
for each one of the others separately. In other words if any three sta- 
tion-position points are col linear, the configuration is singular. To 
derive the determinantal loci of any case formally, one should consider 
the expression of the determinant of a typical submatrix analogous to 
(107), and equate each factor of the typical term to zero. Consider for 
example the case of six stations, ( x i >y± ^0 > (i=l,2,3,4,S,6) observ- 
ing continuously range difference in a pattern of observation such that 
three station quadruples, say 056=(1, 2,3, 4) , 046^(1,2,3,5) and 
Q45 = (l, 2,3,6) carry out all the observations. Let {3 5^ } , {j 45), (j 45} be 
the satellite-position point sets observed by 056, 046, 045 respectively. 
Then the determinant of a typical submatrix of the condensed jacobian 
matrix of this case may be written as follows. 
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! A (°S6> 0 46’°45) I" 

P 1TP2 (p 3 


(-l)°l +0 2 +a 3 +a 4 a 10 


(P4 (p 5) ) ) ) 


b 2m2 1 


a 2m2 D 2m2 

a 2m3 b 2m3 


x o^4 » m 5 » w 6d™ 7 * m 8 > m 9p ™1 0 » m l 1 > m l 2 d™ 1 3 » m l 4 » m l 5 

mi m2 m3 1114 ms 1% 
x A la6 A 2y6 A 2e^ A 3n0 A 3KX A 3v5 

where 

m 1 ,m2,m3,m4,m 5 ,m 6 e{j56} U { 145 ) U { J45 ) 
m 7> m 8, m 9 £ ^ 56^ 
m 10> m ll» m 12 e ^ 46) 

< n 13» m 14» m 15 e ^ 4S ^ 

1 (56) ,if mj ,m2,m3 ,ni4 ,ms ,mge { jgg ) 
(46), if m 1 ,m 2 ,m 3 ,m4,m 5 ,m 6 e{j4 6 } 
(45), if m 1 ,jn2,m 3 ,m 4 ,'m5,m 6 e{j45}- 


(116) 


(117) 


a so 

b so 

c so 

1 



a smi 

b smi 

c smi 

1 

, (s=3,4,5,6) 

(118) 

a smj ( 

bsmk 

Csm k 

1 



a sm£ 

b smj, 

c smj. 

1 




The determinantal loci are obtained by searching for the zeroes of the 
typical term of the summation (116). 

Determinantal loci for the satellite-position points . 

(i) From 

a 10 1 

= 0, (m 1 e{j 56 } U{j. 46 ) U{j 4 s>) 
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it is concluded that the circular cone with vertex at PjE ( xj^ , , axis 
parallel to the x-axis and passing through the point Q c = (X Q , Y Q , Z Q ) is a 
determinantal locus for the satellite-position points (jg^ } U{j 46 J U{ j 45 } • 

(ii) From 

a 20 b 20 1 

a 2 m 2 b 2 m 2 1 “ °> ( m 2 > m 3 e < 156 * U(j 46 > UtlUstt 

a 2m3 b> 2m3 * 

it is concluded that any two semi-straight lines through P 2 = (x 2 ,y 2 » z 2 ) » 
(which may coincide) constitute a determinantal locus for the satellite- 
position points (i 5^} uf j 4 gJ U{ ^45} • Since the two semi -straight lines 
may coincide, any single semi-straight line through P2, is also a deter- 
minantal locus for the satellite-position points {j 46 ) {j 45 }. 

(iii) From 
nu.mc: .nu; 

d 3 * 0, (m 4 ,m 5 ,m 6 e{j 56 }u{j 46 )U{j 45 } ) 

it is concluded that any three semi-straight lines through P3= (X3 ^3, Z3) , 
(any two of which or all three may coincide), constitute a determinantal 
locus for the satellite-position points (j 55) L){j 4 g} Ut 145^- 

(iv) From 

m 7 j R j Til g 

°4 = °> ( m 7 > m 8 > m 9 e O 56 -^ 

it is concluded as above that, any three semi-straight lines through 
P4=(x 4 ,y 4 , z 4 ) , (any two of which or all three may coincide) constitute a 
determinantal locus for the satellite-position points (15^), 

(v) From 

n mio>mn,m 12 ...... 

d 5 = °» ( m 10 ' m ll> m 12 e{ ! 46 } ) 

it is concluded that, any three semi-straight lines through p 5= ( x $ ,y 5 , Z5), 
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(any two of which or all three may coincide) constitute a determinants! 
locus for the satellite-position points { j 4 6^ - 

(vi) From 

D ™13» m 14» ra 15 = 0, (m 1 3,m 14> m 15 e{ j 4 5 }) 

it is concluded that, any three semi-straight lines through p 6= ( x 6»y6> z 6l. 
(any two of which or all three may coincide) constitute a determinantal 
locus for the satellite-position points { j 4 5 1 . 

(vii) From 

I m l e ^56^ U ^46^ (J ^45^ 

1(56), if m 1 c{j 56 } 

) 

(ag)=((46), if 

((45) , if m 1 e{j 45 } 

it is concluded that, the planes ( p 2> p 3> p 4)» ( P 2 ,P 3 ,P 5^ anc * ( P 2 ,P 3 ,P 6^ 
constitute simultaneously a determinantal locus for the point-sets { j 56 * 
(j 46 J and { j 4 s> respectively. 

(viii) From either one of the relations 


A 2y6= °» A 2ec = °*j(Y«)»(cO = 


m 2 ,m3e{j S 6}U{.i46J U<j45> 

! 

(56), if m 2 ,m 3 e{j 56 } 

(46) , if m 2 ,m 3 e{j 46 } 

(45), if m 2 ,m 3 E{j 45 } 
it is concluded that the planes (Pi,P 3 ,P 4 ), (Pi,P 3 > p 5 ) and ( p i> p 3> p 6) con- 
stitute simultaneously a determinantal locus for the point-sets {jsg}, 
{j 46 }, and { j 45 } respectively. 
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(ix) From either one of the relations 

,n 4 » m 5> m 6 e {^ 56>IJ^46^ U ^45^ 

1 (56), if m^,mg,m^c{ j ^.^respectively 

(46), if m 4 ,m 5 ,m 6 e{j 46 }respectively 

(45), if m^.inj.mgetj^^lrespectively 

it is concluded that the planes (Pi>P 2 >P 4 )> (Pl>P 2 »^ 5 ^ an ^ 

constitute simultaneously a determinantal locus for the point-sets 

(j 46 ) and {j 45 ) respectively. 

(x) From either one of the relations 
my me mg 

A 456 " » A 456 " °» A 456 = °» (>"7,m 8 ,m g e{ j 56 >) 

it is concluded that the plane Pi>P 2> p 3 is a determinantal locus for the 

point -set (js^)- 

(xi) From either one of the relations 

min m ll mi 2 

a 546 3 A S46 = a 546 = °> ( m 10 > m \ 1 > m 12 G 1 J 46^ 

it is concluded that the plane (Pi,P 2 > p 3 ) i s a determinantal locus for the 
point-set {j 46 }. 

(xii) From either one of the relations 
m 13 m 14 mjs 

a 645 = °> a 645 " °> a 645 " °> ^13^14^15^^ 45 » 
it is concluded that the plane ( p i> p 2 , ^ 3 ^ a determinantal locus for the 

point set { j 45 } • These are the determinantal loci for the satellite- 
position points. 

Determinantal loci for the station-position points 

They are obtained by searching for the linear dependence of the last 
three rows of the determinants v in the expression (116). It seems 
easier to consider these determinants back in their positions in the con- 
densed jacobian matrix. That is not given in this case, for it is of the 


m 4 

A 3n<T°’ A 


m 5 

3<X 


n>6 

=°> A 3vr 0 ’ 
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same form with the corresponding one for range networks (see Table 9), 
and the following correspondence holds between their entries. 

(Xj- x i)A^u V ^f a ij “ a i 0 )^inv 

(Yj-yiWpv ' <_4 '(t > ij *t>i 0 ) Axjjv ( - 119 - ) 

(Zj-Zi)Axyv ‘* -> ( c ij _c io^4u v 

Then one has the following. 

(i) By considering the determinants 

4 56* A 146’ A 145’ ( for an Y ^ > 

of the first three columns (xi ,yi , zj-column block) of the condensed jaco- 
bian matrix (refer to Table 9), one concludes that, if each one of the 
point-triples 

{P 2 > P 3 » P 4 >> { p 2> p 3> p s}> ( p 2* p 3> p 6} 

is collinear, the first three columns (xi ,yi , z \ -column block) of the con- 
densed jacobian matrix become zero ones. Therefore the straight line 
determined by ?2 and P 3 is a determinantal locus for p 4, p 5 and p 6 - 

Similarly from the rest three-column blocks of the condensed jaco- 
bian matrix one obtains the following. 

(ii) From the x 2,y2> z 2" coluTnn block it is concluded that, if each one of 
the point-triples 

{ p l,P3, p 4>* { p l, p 3’ p 5>> f p lg P 3.V 

is collinear, the X2 ,Y2> z 2“ column block becomes a zero one. Therefore 
the straight line determined by Pi and P 3 is a determinantal locus for 
P 4 ,P 5 and P 6 . 

(iii) From the xj^, zj-column block, it is concluded that, if each one 
of the point-triples 

{ p l, p 2 » p 4 ^ { p l. p 2» p 5>» { p l» p 2 * p 6 > 
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is collinear, the X3,y3, Z3-column block becomes a zero one. Therefore 
the straight line determined by Pj and ?2 a determinantal locus for 
P4 , P 3 and P 6 . 

(iv) From the last three column blocks, i.e. 
x 4»y4> z^column block, 

X5,y5» z^-column block, 
x 6»> r 6 ,z 6“ colu,1in 

it is concluded that, if the points p i, p 2> p 3 are collinear, all these 
three column blocks become zero ones. Therefore the straight line deter 
mined by any two of the three points p i, p 2» p 3 is a determinantal locus 
for the third one. 

Remark : In this observational pattern, it may happen that the point- 

triples above have two points in common. However this is not the 
case for all the observational patterns. In order to find out 
the general form of a determinantal locus of this kind, one should 
notice that it is associated with one of the stations and is defined 
through the collinearity conditions of the point-triples which result 
from the station quadruples in which the associated station participates 
after excluding that station. The general procedure of finding out such 
a determinantal locus is given in the following section 2 . 4 . Recapi- 
tulating the results of the above observational pattern one has the 
following. 

Determinantal loci for the satellite-position points 
( 1 ) Any three semi-straight lines (any two of which or all three may 
coincide) through either one of the station-position points P 3» P 4* P 5» P 5 
(with no coordinates constrained) constitute a determinantal locus 
for the satellite— position points observed by the station-quadruples, 
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in which the corresponding station participates. 

( 2 ) Any two serai -straight lines (which may coincide) through station P2 
(whose one of the coordinates, z 2 , is constrained) constitute a determm- 
antal locus for the satellite-position points observed by the station- 
quadruples, in which the station participates* 

( 3 ) The circular cone with vertex at the station-position point Pj (when 
two of the coordinates, yj.Zi, are constrained), axis parallel to x-axis 
and passing through the point Q 0 = (X 0 ,Y 0 , Z 0 ) , is a deterrainantal locus for 
the satellite-position points observed by the station-quadruples, ... 

in which the station P participates* 

( 4 ) The plane (P1.P2.P3). is a determinantal locus for each one of the 
point-sets {356}. (j4f,)» O45) separately. 

( 5 ) The planes (P 2 .P 3 .P 4 ). (P 2 .P 3 » P 5 ). (P2. P 3.P6) constitute simultane- 
ously a deterrainantal locus for the point-sets (ish)* ^46^* ^i45^ respec 
tively. 

(6) The planes (P 1 ,P 3 ,P 4 ), (Pj.P3.P5). WlWs) constitute simultane- 
ously a deterrainantal locus for the point-sets (j 56 }> ^ 46 *’ ^45^ res P ec " 
tively. 

( 7 ) The planes (P1.P2.P4)* C P 1 > p 2 > p 5 ) > (Pl.P2.P6) constitute simultane- 

ously a determinantal locus for the point-sets i j 56 ^ » 046 ^* ^ 45 ^ res P ec " 
tively. 

Determinantal loci for the station-position points 

(8) Each one of the straight lines (Pi,P 2 ). ( p 2 > P 3 ) ( P 1 ’ P 3 ) is a deteT " 
minantal locus for P4,Ps and P 6 . 

( 9 ) The straight line determined by any two of the three points P1.P2.P3 
is a determinantal locus for the third one., ... 
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Remark: The selected set of admissible constraints 


X 0 ,Y 0 ,Z 0 ,yj, ,22 - constant appears to affect the determinant a 1 loci in 
(2) and (3). Thus, as in the case of range-networks, the selection of 
the set of admissible constraints does affect the determinantal loci. 

2.4 Synopses 

After the above discussion one concludes the following about the 
determinantal loci for networks of ground stations observing range dif- 
ferences by means of geoceiver. 

Determinantal loci for satellite-position points 

(а) Any three semi-straight lines (any two of which or all three may 
coincide) through a station-position point P^s(x^,y£, z^) , ( with 

no coordinates constrained) constitute a determinantal locus for the 
satellite-position points observed by the station-quadruples, in which 
station P| participates. 

In case one of the coordinates of Pj is constrained, the three semi- 
straight lines are reduced to two, while in the case when two of the coordinates 
of Pi are constrained, the determinantal locus is the surface of the cir- 
cular cone with vertex at Pj , and axis parallel to the coordinate axis 
corresponding to the unconstrained coordinate. 

(б) A second category of determinantal loci for satellite-position points 
consists of sets of planes determined by the station-position points, (see 
determinantal loci (4), (5), (6), (7) at the end of the preceding section 
2.3). The determination of a determinantal locus of this category, de- 
pends upon the observational pattern^ however one can give a general pro- 
cedure of determining a determinantal locus of this category for any 
observational pattern, without having recourse to the condensed jacobian 
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matrix, (see following Remark). Each determinantal locus of this cate- 
gory is associated with a ground station. In other words there are as 
many determinantal loci of this category as ground stations. Consider 
any station of the network, say Pj . The determinantal locus associated 
with this station is determined as follows: 

Let Oi,02,...,Ofc be the station-quadruples in which participates, 
and {Qj } i , )k be the sets of satellite-position points ob- 
served by these station-quadruples respectively. Excluding station- 
position point Pj, the remaining three station-position points in each 
one of the station-quadruples Oj, 0 2 > • • • jO^ determine a plane. Let these 
planes be denoted by tt},it 2 , . ■ . Then the planes irj ,it 2 , . . • consti- 

tute simultaneously a determinantal locus for the satellite-position 
point-sets {Qj ) lf {Qj • ■ - , (Qj \ respectively. That is, if 

(Q j } l e1 Tl 
{Qj>2 eiT 2 

{Qj>k eir k 

there is no unique least squares solution. 

One may check the procedure with the observational pattern of the 
preceding section 2 . 3 . In that observational pattern the three station- 
quadruples 055,045 and O45 carry out all the observations, observing the 
sets of satellite-position points 56> - 1 J 46 > and {145} respectively. 

The determinantal locus of the category under discussion, associated with 
station Pj for example, according to the above described procedure, is 
determined as follows. The station participates in all three station- 
quadruples 055,045 and O45. By excluding , the remaining three station- 
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position points in each one of the station-quadruples 0 56 ,0 46 ,0 45 deter- 
mine the planes 

CP 2 ,P3. P 4^> ( P 2> P 3> P s)« ( P 2> P 3> P 6) 

respectively. These planes constitute simultaneously a determinantal 
locus for the point-sets {j 56 } , { j 46 h U 4S> respectively. This checks with 
the result of the preceding section 2.3. As a second check, take the 
determinantal locus associated with the station P 4 . This station partic- 
ipates in the s tat ion -quadruple 0 56 =(1,2,3,4) only. Bv excluding P 4 , the 
remaining three station-position points determine the plane 

( p l> p 2’ p 3) • 

This plane is a determinantal locus for the satellite-position points 
{j 56 } observed by the station-quadruple 055 = (1 , 2, 3, 4) . This result checks 
also with that of the preceding section 1.3. 

Remark: In case of a determinantal locus of this category, say that 

associated with station P^, the three columns corresponding to the coor- 
dinates of Pj, of the condensed jacobian matrix become zero ones. 

Determinantal loci for station-position points 

Each of these determinantal loci is associated with a station. In 
other words there are as many determinantal loci of this type as stations. 
The general procedure of finding out such a determinantal locus is analo- 
gous to that which was described above for the satellite-position points. 
Consider any station of the network, say P^. The determinantal locus for 
station-position points associated with this station is determined as 
follows. 

Let 0^,02 ,..., 0^ be the station-quadruples in which participates. 
By excluding station-position point Pj_, there remains from each one of 
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the station-quadruples 0 1 ,0 2> . 
be denoted* by 

P|,P2.P3 ■ P l- P l> P i ■ ' 


. ,0 k a point-triple. 



P 



k 

3 • 


Let these triples 


Let also 

(cj) » C e 2^ * • * • » ( e k-* 

be the straight lines determined by any two of the three points of the 
above point-triples respectively. Without loss of generality one may 
take 

( ei ) = (pJ,P2), (e 2 D=( p l.? 2 )» (ek) = ( p l» p 2) 

Then the straight lines (ej) ,(e 2 ) , . . . , O k ) constitute simultaneously a 

12 k 

determinantal locus for the station-position points P 3 »• • • > p 3 respect- 
ively. That is, if each one of the point-triples 

mS.pJ.pJ). {p?.p1.pIj. .... cp5.p5.pSi ' 

is col linear, the Xj^.y^.z^-column blockt of the condensed jacobian matrix 
becomes a zero one, and there is no unique least squares solution. 

Remark: It is concluded from relation (97) that, if one knows the n 
ranges S? (i=l, 2, . . . ,n; n=the number of the stations), then the problem is 

reduced to one of ranges. 


*For a specific observational pattern these triples are completely deter- 
mined. The above notation is used in order to facilitate the general 
discussion for any observational pattern. 

txi,yi,Zi are the coordinates of the station-position point Pi. 
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APPENDIX A 

Implicit Function Theorem 

The implicit function theorem is stated here without proof. It is 
given in the same form as in [Caratheodary, 1965], (p.10). 

Let n functions f 1} f 2 ,f n be given depending on the variables 
Xj (j = l,2,. . . ,n) and t 0l (a=l ,2 . . ,m) and being continuous together with 

all their first partial derivatives 


8f-i , 9fj ( A1 ) 

3xj 3tot 

in the (n+m) variables (Xj,t a ). 

At an interior point (Xj, t a ) in the domain of definition of f-^ the 
equations 

fiCXj, t*) = 0, i=l,2,...,n (A2) 

are all satisfied and also, the determinant 


O 



(A3) 


There then exists in a certain neighborhood of the point (t a ) one and only 
one system of continuous functions 


Xj= 4>i Ct a ) , (i=l ,2 . ,n) , (A4) 

which, first, satisfies the conditions 

X- = ^(O, (i=l,2,...,n) (A5) 

and second, satisfies identically the relations 

. f i (<fri (t qj) , t g) = 0, (i, j=l, 2 , . . .n; , =1 , 2 , . . . , m) . (A6) 

The first derivatives 94 i ^/8t a exist in the same region, are continu- 
ous functions of t a and are found by solving the system of equations 
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(A7) 


Mi Mi + Mi = o, (i=i, 2 ,...,n). 

9Xj 9ta 3 tot 

As an example consider the functions of the left-hand sides of ( 1 ), for 
i-1,2,3. 

f i (Xj ,Yj , Zj ,xj .yj » 2 j) = [ (Xj-x 1 ) 2 +CYj-y , i) 2 + (Zj-zi) 2 ] 
f 2 (X j ,Y j ,Z j ,x 2 ,y 2 ,z 2 ) - [(X r x 2 ) 2 + (Y r y 2 ) 2 +(Z r z 2 ) 2 ] (AS) 

f 3 ( X j* Y j’ 2 j * x 3 ^ 3 > z 3 ) = [( x j-x 3 ) 2 +(Y ; j-y 3 ) 2 +(Zj-z 3 ) 2 ] 

The correspondence in notation between this example and the statement of 
the theorem above is as follows. 

x j, Yj, Zj, Xj , y 1} zj, x 2> y 2 , z 2 , x 3 , y 3> z 3 

X l» X 2» X 3» t l» t 2» t 3 J t 4» t 5> t 6 * t 7’ t 8 ’ t 9‘ 

The above functions are continuous with continuous first derivatives; in 

fact they are analytic. Let P= (Xj , Yj , Zj , xj ,y| , , x 2 , y 2 , z 2 ,x 3 , y 3 , z 3 ) be a 

point which annuls f 2 ,f 2 , and f 3 , i.e., 

fj (Xj >Yj , Zj >XjVYi > z P = 0 

f 2^ X j » Y j » z j » x 2 » y|» z |) = 0 (A9) 

f 3 CX],Y?,Z?,x^y^,zp * 0 
and suppose that 

x r x i Y j-yi Z r z i 

x r x 2 y]-y°2 z r z 2 (aio) 

Xj-x§ Yj-y| ZJ-*S 

Then according to the above theorem in a certain neighborhood of the point 
( X 1 z \> x 2 x 3 »^ 3 ' Z 3 ) there exists one and only one system of con- 
t inuous functions' 
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(All) 


Xj = <l>i(x 1 ,y 1 ,z 1 ,x 2> y 2 ,z 2 ,X3,y 3 ,Z3) 
Yj = ^(xj, |,Zi>x 2 ,y 2 ,z 2 ,x 3 ,y 3 ,z 3 ) 
Zj = 4> 3 Cx! J y 1 ,z 1 ,x 2 ,y 2 ,z 2 ,x 3 ,y 3 ,z 3 ) 
which satisfy identically the relations 

C 1 » 4»2 >^3» x l ‘Xl » z l^ = 0 

f 2 Ui’*2’ ,{> 3 jX 2’> r 2’ z 2) = 0 

f 3^1 >((, 2^3’ x 3»> r 3» z 3^ = 0 
In other words one may solve the system (1). 
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APPENDIX B 

The Intersection of Three Spherical Surfaces 
It is known that two spherical surfaces intersect in a circle* 
the plane of which is perpendicular to the line connecting the center 
of the spheres* Also, the intersection of a plane and a sphere is 
a circle with its center at the foot point of the perpendicular 
from the center of the sphere to the plane* The procedure to be followed 
is to consider the plane of the intersection of two of the three spheres* 
the intersection of that plane with the third sphere, and then on 
that plane the intersection of the two resulted circles* The points 



Section of the three spheres with the plane of the centers 

Figure Bl * 
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intersection (if they exist) of these two circles are the points of inter- 
section of the three spherical surfaces. Since two circles intersect at 
most at two points (real), the set of intersection points of three spheri- 
cal surfaces contains at most two points which in fact, are symmetric with 
respect to the plane of the centers of the three spheres. To formulate 
the above ideas analytically, let (C^), (C2), (C3) [see Figure Bl] be the 
three spheres with centers the ground-station points Pi»P 2» P 3 respective- 
ly. Consider the intersection of (Cj) and (C2) and let t * ie 

circle of intersection, (A its center and pi 2 its radius) and it the plane of 
that circle. Let also C 37r (B,p3 7r ) be the circle of intersection of the plane 
n and the sphere (C3) . For convenience a table of the notations follows. 

Ci (Pf , r i j ) , (i=l,2,3) - Sphere (Cj ) with center P^ and 

radius r^j . 

Ci 2 (A>P 12) " Circle °f intersection of the spheres 

(C.) and (C2) with center at A and 
radius equal to p-^- 

•nr - Plane of intersection of the spheres 

(Ci) and (C 2 ). 

C3i T (B,p 3ir ) - Circle of intersection of the piano 

tt and the circle (C3). 

fi , r 2 > r^j^A* ^B ~ Position vectors of the points 

P 1> P 2» P 3> A » B respectively. 


The analytic expressions of the spherical surfaces are 


(Ci) : 

11*1 

1 2 

2ri *X+ | 

I'll 

|2 _ 2 
1 r l j » 

(Bl) j 

(C 2 ) : 

11*1 

1 2 - 

2r 2 * X+ | 

|r 2 l 

I - r 2j . 

(Bl ) 2 

(C 3 ) : 

11*1 

1 2 - 

2r 3 -X+ | 

1^3 1 

m2 _ _ 2 
1 “ r 3j » 

(Bl ) 3 

where vector 

notation 

is used 

for 

convenience. 
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Notice that 


X - (x,y, z) , r j - ( x i > yi » 2 i ) * ^ _ ^ ? 2 , 3 , 


(B2) 


j | X j 1 2 = (x 2 +y 2 +z 2 ), | | r ± | 1 2 = (x?+yf+z?), i-1.2,3 
and r i *X = (x^+yiV+z^) 

The analytic expression of the plane tt (intersection of the spheres (C x ) 
and (C 2 )) results by subtracting (Bl)-[ from (Bl) 2 , i.e.. 


tt : 2(r 2 -r 1 )X = ( I |r 2 | | 2 ' r 2 j) ~ < I l f 1 1 1 2 ‘ r lj > 


or 


it : t • X . = p , 


where 


p « (II^MMp - (M f ill 2 - r ip- 


(B3) 

(B4) 

(B5) 

(B6) 


In order to find out the position vectors of the centers and the radii of 
the circles C 12 (A,p 12 ) and , one considers the intersection of 

the plane tt with (C 3 ) and (Cj) (or C 2 ) respectively. Since BP^/P^ one 
may write for the equation of the line BP 3 , (see Figure Bl). 

BP 3 E &2 : X=r 3 +A(f 2 -ri]=r 3 + Xt, (X = a scalar). (B7) 

By considering B as the intersection of tt and BP 3 =£ 2 , its position vector 
results from (B4 and (B7) , namely 


X-t = p 
X=r 3 +Xt 


f B -t = p 


r B“^3 +A B t 


Xb = 


P=^ 3 -t 


BP- 


(B8) 


and 


r P-r 3 *t.,_ 
r B“ r 3 + t '. ■ ■ ~ ] t 


(B9) 


lltll 
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For the radius p 3lT one observes (see Figure Bl) that 

1 2 


2 2 
p 3tt - r 3j 


BP 


31 


and by (B 8 ) 

2 _2 rlZllli 2 

P3tt - r 3j E || t j|] • 


P-r i *t 


(BIO) 


Analogously for the circle Cj 2 (A,pi 2 ) one obtains 

r P - ? l* f ,= 


f A =f l + E 


I it 


rlt. 


P-r, • t 


2 r 2 r 1 1 

12 ij " [ 7WT 


(Bll) 


(12) 


by considering the intersection of it with (Cj), or the equivalent 

- - r P_ t2*t n 

r A =r 2 ♦t '■ - — l 1 . 

IMP 

7 7 P — T 9 4 t 7 

''TNI 1 ' 


(B13) 


(B14) 


by considering the intersection of tt with(C 2 ). This is all that is needed 
for the conditions of intersection of the three spherical surfaces. Ob- 
viously the common points to the three spherical surfaces are the common 
points to the two circles, C^fA.p^) and both of which lie on 

the plane it. The conditions for two circles in a plane to intersect are 
well known, i.e., the distance of the centers is greater than the abso- 
lute value of the difference of the radii and less than the sum of the 
radii. Thus the required conditions are 

l p 12“ p 3-J < M f A' f Bil < p 12 +p 3 tt * (B15) 

For one obtains from (B9), (Bll), and (B5) 

i (f r f 3 )-(f 2 -f 1 ) 




l _i 3 2 


l r 2 -r l 


-(r 2 -fj) 


(B16) 
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or from (B9) , (B13), and (B5) the equivalent expression 


= (^2^3) - ^ y- —y (r 2 -ri) 


(B17) 


For (p 12 -P 3 t!) and fPl2 + P3ir ) one uses the formulas CB12), (BIO) or (B13), 
(BIO). It is understood that one may consider first, instead of the 
intersection of the spheres (Cj) and (C 2 ) that of the spheres (C 2 D and 
(C 3 ) or (C 3 ) and (C x ). One observes also that the intersection of the 
three spheres in pairs is a necessary condition but not sufficient. That 


is, the conditions 


1 r l j ~ r 2j ! 

< 1 

l?2’ ? lll < ( r l j + r 2j ^ 

(B18)j 

I r 3r r 2j | 

< 1 

l f 3- f 2M < ( r 2j +r 3j) 

(B18) 2 

l r 3j“ r ljl 

!•< I 

m f ii 1 < ( r 3i +r ij) 

(B18) 3 

are necessary 

for 

(B15) but not sufficient. 

After the conditions for the 


intersection, one may proceed to find the degenerate cases. 

The left-hand side of (3) is the determinant of the matrix whose 

I 

rows are the vectors (fj-rj), (?j-r 2 )> (rj-? 3 ). Thus condition (3) may 
be interpreted as the condition for the vectors (r-j-fj) , (fj-f 2 ) , 
to be linearly independent, or equivalently n£i to be coplanar. This is 
equivalent to saying that the four points Qj,Pi»P2» p 3 are not coplanar. 

Since Qj , as it was said above, does not lie in the plane of the centers 
p l* p 2» p 3» the only case for the points Qi»Pi>P2» p 3 to coplanar is when 
P 1» P 2» P 3 are collinear. Thus (3) is equivalent to stating that the cen- 
ters of the three spheres are not col linear. In the case where the centers are 
collinear, the three spherical surfaces which have one point in common, 
have in common the circumference of the circle through that point whose 
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plane is perpendicular to the line of the centers, and the center lies on 
that line. Thus under conditions (3) and (B15) the set of points which 
annul the system (2) contains two points symmetric with respect to the 
plane of the points Pi»P2** ) 3- 
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APPENDIX C 


Th e Invariance of the Rank of the Jacobian Matrix 
The invariance of the rank of the jacobian matrix is equivalent to 
the invariance of the jacobian determinant. To see that, consider a jaco- 
bian matrix of rank y; that means that there exists at least one yxy sub- 
matrix whose determinant is non-zero; but that determinant is the jacobian 
determinant of the corresponding functions; therefore the rank of the 
jacobian matrix is invariant under some transformations if and only if 
the jacobian determinant is invariant under those transformations. 

Consider the set of equations 

y l = ( C1 ) 

* f n( x l»'--» x n>- ■ 

which may be looked upon either as a transformation of coordinates or as 
a transformation of the points of the vector space v n of coordinates Xj 
among themselves. Let be a point which corresponds to 

(y° » y 2» * * ’ ,y P‘ T ^ en in tlie neighborhood (xj,...,x“), consider the 
differential space B x which is spanned by the differentials dx 1 ,...,dx n , 
and in the neighborhood of (yi»...,yn) the differential space By which is 
spanned by dy 1 ,...,dy n . The jacobian matrix is the matrix of the linear 
transformation under which the differential space B x goes into the dif- 
ferential space By, i.e., 

(C2) 
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or 


dy = J M *dx. (C3) 

Suppose that dxj, dx2,...,dx n are independent from each other, that is 
they constitute a basis for the differential space B x , and this is the 
case for dyi,...,dyn in By. 

Consider the affine transformation 

Z = Ax+B, [ A 1 5^0, (C4) 

where A is a constant nxn matrix, 3 is a constant n-vector, x= (x^ , . . . ,x n )^ 
and z- (zj . . , z n ) T the vector onto which x goes under the transformation 
(C4). Then the jacobian matrix in (C3) is transformed as follows: 

3f^ n 9fi 3zi, n 9f; 

3x]~ = kll 3zk 3xj~ = kll 3zk a kj fC5) 

where a^j is the (k,j) element of the matrix A. Then one may write in 
matrix notation 


J M (x) » J M (z)A. (C6) 

By taking the determinants of both sides one has 

|J M (X)| = [J M (z)|*N = I a M j m( z ) I ■ (C7) 

Hence |j|^(x) | is a relative invariant of weight* one, under the transfor- 
mation (C4) . Differentiation of (C4) yields 

dz = Adx or dx = A~*dz, (C8) 

and then substitution of (C6) and (C8) into (C3) 

dy = J^(x)dx = J m (z)A A'^dz = J M (z)dz. (C9) 


The last result shows that J^(x)dx is an invariant, under the transforma- 
tion (C4) . 

*It is recalled that the weight is rendered by the exponent of the deter- 
minant of the matrix of the linear operator. If the exponent is zero 
then one has an absolute invariant, otherwise a relative one. 
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Since the jacobian determinant is invariant under the affine trans- 


formations, the rank of the jacobian matrix is invariant also, under the 
same transformations. 
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APPENDIX D 

Pivotal Condensation ([5], pp* 45-48) 

By condensation of a determinant it is meant the reduction of its 
order in a systematic way. There are various methods of condensation, 
from which the method of pivotal condensation is presented here. Since 
pivotal condensation corresponds exactly to the process of successive 
elimination which is used in the solution of simultaneous linear equa- 
tions, it seems to be easier to introduce the method from that side. 

Consider for example the following set of simultaneous equations: 

a l x l + ^l x 2 +c l x 3 +c ^l x 4 “ ^1 

a 2 x l +b 2 x 2 +c 2 x 3 +d 2 x 4 = h 2 (Dl) 

a3Xi+b3X2 + C3X3+d3X4 - I 13 

a4X 1 +b4X2+C4X3+d4X4 = h 4 , 

and try to eliminate the variable Xj. To do that one may proceed as 
follows: by multiplying all the equations except the first by a^ and then 
performing the operations 


2 nd eqt . 

- a^l st eqt) 



3= eqt. 

- a^l st eqt) 



4= b eqt. 

- a 4 0 st eqt) 



to obtain 




a l x l + 

bix 2 + 

01x3+ 

CL 

i— > 

X 

II 

_ cr 


0 + (a 1 b 2 -b 1 a 2 )x2+(a 1 C2-c 1 a2)x3+ (a 1 d 2 -d 1 a 2 )x4 = h 2 -a 2 h 1 

0 + (a ^ b 3 — b j^a 3 ) x 2 + (a 1 c 3— c 1 a 3 ) x (a 1 — d ^a = h^— a^h-^ (D2) 

0 + f a 1^4"bl a 4) x 2 + ( a l c 4~ c l a 4^ x 3 + (a jd^-dja^)x^ = h^-a^h^ 


preceding page blank 
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Now consider the 

determinant 

of the 

coefficient 

matrix , 

and observe 

the 

above operations 

for the 

elimination 

of 

x l> 

on 

that determinant. 





a i 

c 

1 

d l 








ai 

b? c ? 

d 2 





1 A | = 

! a l b 2 

c 3 d 4 | = 

Ls 



, ( 

*) , 


(D3) 



a 3 

b 3 c 

3 

d 3 








a 4 

b 4 c 

4 

d 4 






a l 

b l 

C 1 

d l 



a l 

b l 

C 1 

d l 

a l 1 A | « 

a l a 2 

a l b 2 ! 

a l c 2 

a l d 2 



0 

l a l b 2 1 

i a l c 2 1 1 

a l d 2 1 

a l a 3 

a l b 3 ; 

a l c 3 

a l d 3 



0 

l a l b 3 1 

i a l c 3 1 1 

a l d 3 I 


a l a 4 

a ! b 4 

a l c 4 

a l d 4 



0 

1 a l b 4 1 

l a l c 4 1 1 

a l d 4 | 


|a I = ^ 

a i 


l a l b 2 I 
I a l b 3 I 
I a l b 4 I 


I a l c 2 1 
l a l c 3 I 
I a l c 4 I 


l a l d 2 I 
l a l d 3 I 
l a l d 4 I 


(D4) 


Thus the original determinant |Aj has been reduced to one of third order. 

, n-2 . 

In general a determinant of order n is equal to l/a^ times a determin- 
ant of order n-1, the elements of which are the minors of | A | of order 2 
having as leading element placed in proper priority of rows and col- 
umns. Thus the minor belonging to the 1" and (i+1) rows, 1' and 
(j+l) st columns of | A | will be the element of the condensed deter- 

minant. The element aj is called pivotal element or pivot. The pivot 
need not be the leading element, which may be zero. Any convenient non- 
zero element may be taken as pivot, for one could make any element a 

*The representation of a determinant as a row of its diagonal elements is 
very convenient, in this sort of manipulations, and it is adopted in the 
present appendix. 
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leading element by permuting rows and columns. The pivots which are lead 
ing elements are more convenient. Next to the leading pivots in conven- 
ience come the diagonal elements. 

Not only single elements may be used as pivots, but any non-zero 


minor. Principal minors are preferable. The procedure is the same as 
that of a single element pivot. Thus a determinant j A [ of order n which 
is to be condensed by using as pivot a principal minor of order m (m<n) , 


is equal to (A m ) times a determinant of order (n-m) the elements 

of which are the minors of | A ] of order (m+1) having A m as leading prin- 
cipal minor, placed in proper priority of rows and columns. 

Before closing this appendix, a series of important identities are 

given. 


a ll a 22 a 33 I " a ll 


ana22 


a l l a 32 I 


I a 1 l a 23 I 
i a l l a 33 I 


(D5) 




1 a l l a 22 

! i a ll a 23 1 

l a U a 24 1 

l a l l a 22 a 33 a 44 1 

-1 
= a ll 

1 a l l a 32 

1 l a l l a 33 1 

l a ll a 34 ! 



1 a 1 l a 42 

1 \ a i i a 43 1 

I a ll a 44 1 

l a l l a 22 a 33 a 44 1 



l a 1 l a 22 a 33 1 

l a ll a 22 a 34 

= i a ll 

a 22 1 ' 

1 a l l a 22 a 43 1 

1 a l l a 22 a 44 



(D6) 


(D7) 


a ll a 22 a 33 a 44 a 5S 


a l l a 22 a 33 I “ 1 


I a l l a 22 a 33 a 44 I 
! a ll a 22 a 33 a 54 I 


a l l a 22 a 33 a 44 I 


a ll a 22 a 33 a 55 I 


(D8) 
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1*11- 

” a mm a m*l ( m*l 1 

l a ll- 

,m<-2 1 ■ 

. |a u . 

,a mra a in+l,nl 


1*11- 

‘ a mnPm+ 2^+1 1 

i a 1 1 * 

■ a mm a ra+2,m+2 1 * 

• l a n- 

* a tjm a m+ 2 , n 1 


l»ll- 

’ • a mm a ri,m*l 1 

Ian- 

* ■ a mni a n,n]+2 ^ " * 

.u n . 

* 1 


- 134 - 



APPENDIX n 


Condensation of the Jacobian Matrix 
For Four Stations in Range Networks 
i) Pivotal Condensation 

Let P^Xj^y^, Zj) , P2= (*2 »Y2> z l) * p 3 5 ( x 3*y3 * z 3^ > P 4 = ( x 4 ^4* z 4^ be 
four simultaneously observing ground stations, and Qj = (Xj , Yj ,Zj ) , 

(j=l , 2, . . . ,m) be a set of satellite-position points observed by those 
stations. Then the functional relationships on which the statistical 
model is based, are 

f ij*tXj- x i) 2 *( Y r> , i ) 2 +(Z r z i )2 = r u 

i = 1,2, 3, 4 ( E1 ) 

j = 1,2,... , m ,m>6 

(see Sections 1.2 and 1.3). 

This system contains 

v = 4m>24 equations, ( E ^) 

and 

3m+3x4 = 3(m+4) unknowns, ( E3 ) 

among which six are to be constrained, so that there are 

y = 3(m+2) actual unknowns. (E4) 

Nevertheless while forming the jacobian matrix, (J M ) for the system (El), 
no unknown variable is considered constrained. Thus the jacobian matrix 
of the system (El) is of dimensions 4mx(3m+4). Since (El) is cyclicly 
symmetric with respect to Xj.Yj.Zj and x i ,y i ,z i , while writing the explic- 
it form of the jacobian matrix only the columns for Xj and x i will be 
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written down for simplicity. The columns for YjZj 


understood to 
Thus, one has 


follow immediately those of 
for the jacobian matrix: 


and y ± x ± will be 
and respectively. 


(i.j) 

*1 

... 

x i 

... 

X m 

*1 

*2 

x 3 

x 4 

d.l) 

2(X|-Xj) 


■ 


0 

*2(X 1 -Xj) 

■i — - 

0 

0 

0 


... 

| 

■a 

B 

... 

... 


... 

... 

l.J 

0 

■ 

D 



0 

-2(Xj-xj) 

0 

0 

0 

... 

Blfpi 



■ 

... 

... 

. . . 



(l.B) 

B 


0 

. ..0 

2 

- 2 (V*i) 

0 

0 


(2.1) 

2(X|-X2) 

0. .. 

0 

. . . 

0 

0 

-2(Xi-x 2 ) 

0 


• V 

... 

m 


n 

. * . 

. . . 

. * • 

SSBitB 

B m 

(2,j) 

0 

■ 




0 

0 

-2(Xj-x 2 ) 




. . . 

. . . 

. . . 

. . . 

. . . 

. . * 

. . . 


. « » 

(2,m) 

0 


0 

• • 

.0 

2(X m -x 2 ) 

0 

-2(X m -x 2 ) 

0 

0 

(3,1) 

2(Xi-x 3 ) 



■ 

0 

0 

0 

*2(Xj-x 3 ) 

0 

fl! 




B 


• . . 

. . . 

* . . 

... 

(3J) 

0 

...0 



0 

0 

0 

-2{X r x 3 ) 

0 

... 

... 

... 

... 

B 


... 

. . . 


... 

(3,m) 



0 

•- 

.0 

2(X m -* 3 ) 

0 

0 


0 




0 

| 

0 



0 

-2(Xi-X4) 

B 

... 

. . . 



... 

. . . 


. • * 


(4,j) 

0 

...0 

2(Xj-x 4 ) 

0.. . 

0 ! 

i 

0 

0 

0 

-2{Xj -x 4 ) 

... 

... 


... 

... 

T 

... 


. . . 

• * - 

(4,rn) 

0 

... 

0 

.. .0 

2(V*4)| 

0 

0 

0 

-2 (X m -X4) 


Table El . The Jacobian Matrix for the Case of Four Stations 
Note : Only the x-columns appear for brevity. 
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According to what was said in Section 1.3, the least squares solu- 


tion is unique iff 

rank = y = 3(m+2) 


(ES) 


while for 


rank < p = 3(m+2) 

one has the determinanta 1 loci of the configuration. Therefore in order 
to find out the determinanta 1 loci, one has to consider the determinants 
of all yxy submatrices of J^, where the value of p = 3(m+2) is determined 
through the value of m, which is given from the equation 
rm = 3(m+2) , i.e., m=6 . 


Thus 

p = 3(m+2) = 24 . ( E7 ^ 

It is reminded that m=6 is the minimum number of events for a unique sol- 
ution in the case of four stations, (see Section 1.1.3 relations [6]). 

At this point one is compelled to constrain six of the unknown vari- 
ables, in order to be able to pick up a 24x24 submatrix. However it is 
desirable to arrive at general expressions, before one specifies the 
coordinates to be constrained. In the following the method of pivotal 
condensation will be used to simplify the determinants of those pxp=24x24 
submatrices. Working out that method with a specific set of admissible 
constraints a result will be obtained the symmetry of which will lead to a 
general expression. 

Assume that xj ,yi, zi,y2»Z2> z 3 are constrained and consider, without 
loss of generality, the submatrix of the p(=24) first rows of the jacobian 
matrix. In order to get this submatrix, put m=6 and ignore columns X} ,yi , 
z \*Y2> z 2> z 3 i n Table 1. of the jacobian matrix. Thus the following 
determinant results, where the factor 2 has been dropped from each entry: 
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Table E2 . Determinant of a 24x24 Submatrix of the 
Jacobian Matrix, (see Table El). 
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This determinant should be expanded and its zeroes be searched for the 
ilctcrmi.na.ntal loci. There are many methods of expanding a determinant, 
the most appropriate one here being that of pivotal condensation, (see 
Appendix D) . 

Take the principal minor of the first eighteen rows as the leading 
pivot (assumed non-zero) and let it be called |Ajg[. Also interpret the 
symbol . 

I ^i g , (p jA ) | 

as the principal minor which results by juxtaposing the pth-row and 
A^-column with the principal minor j A 1S | . Then one has for the above 
determinant |&j: 


1 * 1 - 




|Ai 8 , 09.19} |A l8 , {19 p 20)j, 

Kg. (20.19)1, Kg, (2O.20)|. 
Ui8.(21.19)|, Kg, (21 ,20) I , 
(4,8.(22,19)1. 1*18. (22.20)), 
1*18.(23,19)1, |A 18 ,C23,20)|, 


1*18. (2*. 19) | , |4jg, (24, 20) I . 


1*18.(19, 21)|. |A, b ,( 19,22>|, 1 4 lg , (19 .23) | , J A,g , (19, 24) ( 

1*1$, (20, 21)|, |A lg , (20,22) | . Kg- (20.23) I. I* 18 . (20,24) | 

1*18.(21,21)1, [418,(21,22)1, (4,8,(21.23)1. |a 18 , (21,24) | 
1*18.(22.21)1, |4jg, (22,22)), [4,8, (22, 23) j , Kg, (22 ,24) | 
1*18.(23.21)!, |Aj 8.(23,22)1. 14,8,(23,23)1. |4,g, (23,24) ( 
1*18.(24.21)1, (4,8.(24,22)1, |4,g, (24,23) | , (4, 8 , (24 ,24) | 


(E8) 
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This is a compound determinant (i.e., a determinant having determinant 

elements) . For its explicit form the determinant elements are 

expanded by Laplacian expansion: 

(a) Laplacian expansion of [ A ^ 8 I • 

Expanding |A^g[ according to its first three columns there is only- 

one non-zero 3x3-minor which being multiplied by its cofactor gives the 

only term of the expansion. The same is repeated with the cofactor and 

the process continues similarly, until one obtains: 

6 

I Ais! = (-1)°A4-A4 -aJ-A^a5-aJ = -Hj aJ, (E9) 

where 


a=5(l+ 2+3) + (1+7+ 13) + (1+6+11) + (1 + 5+9)+ (1+4+7)-*- (1 + 3+5) = 105 


and 


4 = 


Xj “ X 1 

yj -yi 

Zj -Z 1 


x i 

Y j 

z i 

1 


X 1 

y i 

Z 1 

1 

X r x 2 

Yj ~y 2 

Z j ~ z 2 


x 2 

Y2 

z 2 

1 

Xj-x 3 

Yj ~y$ 

Z j " 2 3 


x 3 

y3 

2 3 

1 


(E10) 


It is appropriate at this place to introduce the determinant 


A^ = 


tj = 


r x i 

Yj-yi 

z i 

- z i 

1 

r x 2 

Y r y 2 

z j 

-Z 2 

1 

j -x 3 

Yj ~y 3 

z j 

" z 3 

1 

X 

1 

Y j -y 4 

z i 

-z 4 

1 

d column manipulations 

A^ 

x i 

yi z i 

i 



x 2 

y 2 z 2 

i 



x 3 

Y3 z 3 

i 



x 4 

y 4 z 4 

i 




(Ell) 


(E12) 
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3 

and A 4 is given by (E10). 

(b) Laplacian expansion of a determinant -element of |a|, say |Ajg, (19, 19) | 
Expanding this determinant according to the last column, only one 
term is non-zero, namely that which corresponds to -{X}-X 2 ). The other 
five terms are obviously zero. To see that, take the minor of -(X 2 -X 2 ), 
and try to obtain its Laplacian expansion according to the columns 4,5,6; 
it is zero, for there is not any non-zero 3x3-minor there. The minor of 
-(Xj-X 2 ) is of the same form with |Ajg|. As a matter of fact, it results 
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from | A 18 | within sign by deleting the 7-th row and introducing the 19-th 
one; in this case the sign changes for the 19-th row get the place of 
the 7-th row by 11 row-interchanges. Therefore working as before one ob- 
tains : 

|a 18 , (19,19)1= - (X 1 -x 2 )&2-A4-A4*aJ-A4-A4 (E15) 

or 

1 6 1 

|a 18> (19,19)1= -(Xj-x 2 ) S .n = + (Xj-x 2 )-Zi |a 19 [ . (E16) 

A 4 A 4 

In an exactly similar way one obtains the expansions of all the de- 
terminant-elements of | A | . Dividing each row of the determinant | A [ by 
| A lg | , and taking out of each row j the common factor I/A4 one has 

-(Xj-x 2 )A 2> (Xj-X 3)A3, (Yj-y3)A3, -(Xi-X4)A 4 , -(Yj-y4)A4, 

-(X 2 - x 2 )A 2 , (X 2 -x 3) A3, CY2-y3)A3, -( x 2" x 4) a 4* -(Y 2 -y4)A4, -(Z 2 -z 4 )A 4 

- (X 3 - x 2 ) A 2 , ( x 3- x 3)A 3 , (Y 3-y 3 )A3 , -(X3-X43A4, -(Yj^^, -( Z 3-2 4 ) A 4 

I A | = (E17) 

4 4 4 4 4 4 

-(X4-x 2 )a 2 » ( x 4 _x 3) a 3» (Y 4 -y3) A3 , -(X4-x 4 )A 4 , -(Y4-y 4 )A 4 , -(Z4-Z4)A 4 

-( x 5 _x 2^ a 2» C x 5 -x 3) a 3* ( y 5"Y3) a 3> ~ ( x S -x 4) a 4 > - ( y 5 _ >'4) a 4> -( z 5 _z 4) a 4 

_ C x 6 _x 2^ a 2> C x 6 -x 3^ a 3> t Y 6~y3^ A 3> '■C X 6 -X 4) A 4> -( y 6'X4) a 4» *( Z 6 _z 4^ A 4 

The symmetry of the above expression suggests the general form one 
is looking for, where the constrained coordinates are not specified. That 
general expression is given below, and it is a sort of condensed jacobian 
matrix . Each time a set of admissible constraints is selected, one has 
to delete the corresponding columns from the general array which follows. 
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Yl 
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nr 

Y2 1 *2 

XJ | ^ 


*4 

*4 

*4 

♦ 

♦ 

<Y r y,)4j 

♦ 

[Zj -Z j ) Aj 

1 

1 j ! 

j (X j“ x 2^2 

r-rt . 1 
1 

^r y 2 ) ^! (Z i- , 2 ,A 2 i 

* ! * 

! ,! 

:tX r *j)6j:(Y r yj)4 J j 

j . !. ... ... 

♦ 

j tZ j" z 3 )A 5 

(Xj-* 4 )A 4 

(Vj-y 4 ^4 

(Z r I 4 )A^ 


j’TT.3. . 


Table E3 . The Condensed Jacobian Matrix for the Case of Four Ground Stations 

For example the expression (E17) of the determinant |a| is obtained 
immediately from the condensed jacobian matrix by deleting the columns of 
the constrained coordinates »yi> z i>y2> z 2> z 3 and selecting for the index 
j the values 1,2, 3, 4, 5, 6 out of the m. 
i i ) By the method of elimination 

It was said in Appendix D that pivotal condensation corresponds 
exactly to the process of successive elimination which is used in the 
solution of simultaneous equations. This fact will be exploited to obtain 
the condensed jacobian matrix by using a technique different from the pre- 
vious one. Consider the differential form of the equations (FI). 

2(X j -x i KdX j -dx i ) + 2(Y r y i )CdY r dy i ) + 2(Z ;j -z i )(dZj-dz i ) = 2r i;j dr i;j . (E18) 

This is a linear mapping of the differential space spanned by the differ- 
entials (dXj ,dYj ,dZj ,dx^ ,dy^,dzj) into the differential space spanned by 
dr^., (i=l,2,3,4 and j =1 , 2 , . . . , m) . Since one is interested here in the 
rank of the matrix of the above linear operator, he may consider the ker- 
nel or null space of that operator, i.e., he may take 

(Xj -xf ) (dXj -dx^ ) + (Y j -y i ) CdY j -dy j ) + ( z,j - zj ) (dZ j -dz — 0, (E19) 

to obtain all the information about the rank. 
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One may group the 4m equations (E19) in m groups of four, of the 

form 

(Xj-XjHdXj-dxjMYj-y^CdYj-dyjMZj-ZiHdZj-dZi) = 0 

fX,-Xo)(dX,-dx 9 ) + (Y;-y 9 )(dY;-dy 7 MZ.-z 7 HdZ.-dz 2 ) = 0 
' J *'■ ■ J J “ ■> ■’ (E20) 

(Xj-xsHdXj-dxsD+CYj-yjHdYj-dysMZj-ZsKdZ^-dzj) = 0 

(Xj -x 4 ) (dXj -dx 4 )+ (Yj -y 4 ) (dYj -dy 4 ) + (Zj -z 4 ) (dZj -dz 4 ) = 0. 

From the above four simultaneous linear homogeneous equations one may 
eliminate the three differentials dXj,dYj,dZj. The resultant of the elim- 
ination is a linear homogenous equation of the form 

a 1 dxi+a2dyi+a3dz 1 +a4dx2+a5dy2 + a.6dz2+a7dx3+a8dy3+agdz3+ai0dx4+aiidy4+ai2dz4=O, 

where the coefficients a^ (i=l, . . . , 12) are the entries of the j-th row of 
the condensed jacobian matrix as it is given in Table E3. The elimination 
of dXi ,dY-j ,dZ.j is carried out formally; one takes three of the equations 
(E20) , solves them for dXj,dYj,dZj and then substitutes in the fourth one to 
obtain the resultant equation. 

This method is followed in (Blaha l 1971]) and (Rinner [1966]) with 

the difference that the constrained coordinates are specified from the 

beginning, and (E20) becomes 

Xj dXj + YjdYj + ZjdZj = 0 

(X-: -x 2 ) (dX^ -dx 2 ) «■ Y.dY; + Z-jdZ.: = 0 
J J J J J J (E21) 

(Xj-x 3 ) (dXj-dx 3 ) + (Y j -y 3 )(dYj-dy 3 ) + ZjdZj = 0 
(Xj-x 4 )(dXj-dx 4 ) + (Yj -y 4 ) (dYj -dy 4 ) + (Zj -z 4 ) (dZj -dz 4 ) = 0 
under the constraints xi=yi=z^=y 2 =z 2 =Z 3 = 0 . 

One observes that symmetry is gone in (E21) , and that fact has many 
consequences for the whole treatment. It is not advisable for a 
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theoretical treatment to destroy symmetry of expressions, for such an act 
deprives one of flexibility in the manipulations. 

The expressions (K14) for the particular set of constraints. 

x^=y 1= z i=y 2 = z 2 = z 3=° > become : 

1 
1 

1 

(E22) 

A 3 = Z j x 2^ 4 ' z 4 x 2y ; j 
A 4 = Zj x 2 y 3> 



x 2 

0 

1 



X i 

Al = Zj 

x 3 

y3 

1 

-z 4 

x 2 


x 4 

^4 

1 



x 3 

i 

x 3 

Y3 



x j 

Y i 

A 2 = Zj 



+■ ! 

z 4 




x 4 

Y4 



x 3 

y3 


^3 
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APPENDIX F 

The Complete In variant System Associated 
With a Set of Geometric Points 

The purpose of this appendix is to quote some applications of the 
theory of invariants to geometry, which are of interest in this work. 

The discussion is within the frame of affine geometry. At first the def- 
initions of fundamental invariants will be given. 

(i) One Dimensional Space 

Consider two points P]^ and P 2 with affine coordinates (x^,l) and 
(x 2 ,l) respectively. Take the determinant 

X 1 1 

A 12 = = xj-X2. CFi ) 

x 2 1 

This two-rowed determinant which is formed from the coordinates of the 
two points is a fundamental invariant in one-dimensional space. Actually 
A|2 an expression of the distance between the two points. 

(ii) Two Dimensional Space 

Consider three points Pi,P2> p 3 with affine coordinates (xj,yi,l), 
(x 2 ,y?,l), respectively. The three-rowed determinant which is 

formed from these coordinates, i.e., 

X 1 Y i 1 

x 2 y 2 1 , CF2) 

x 3 ^3 1 

is a fundamental invariant in two dimensional space. It is well known 
that A^23 is twice the area of the triangle (P1P2P3). 
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a 123 = 


X 1 

Xl 

Z 1 

x 2 

yi 

z 2 

x 3 

y 3 

z 3 

x 4 

> r 4 

z 4 


(iiij Three Dimensional Space 

Consider four points p i » P 2» P 3* P 4 with affine coordinates (Xj /j.Zpl), 
(x 2 ,y 2 , z 2 ,l) (x 3 ,y 3 ,z 3 ,l), (x 4 ,y 4 ,z 4 ,0 respectively. The four-rowed 
determinant which is formed from these coordinates, i.e., 

IXi V i Zi it 


A1234 = 


is a fundamental invariant in three dimensional space. It is known that 
^1234 i s s i x times the volume of the tetrahedron (PiP 2 P 3 P 4 ). 

The theory of invariants proves: 

i) Given a finite number of points on a straight line, 

[Pj=(xj ,1) , . . .P n (x n , 1) ] , the complete invariant system associated with 
these points is the totality of the two-rowed determinants, 
x- 1 i=l, 2, . . . ,n 

Aij - (F4) 

X-| 1 j — l,2,...,n. 

That means that any invariant associated with these points can be ex- 
pressed as a polynomial in the Aij • 

ii) Given a finite set of points in a plane 

[PjrCxj ,yj , 1) , . . . ,P n = (x n ,y n , 1) ] , the complete invariant system associated 
with these points is the totality of the three-rowed determinants. 


x^ y^ 1 i=l,2,...,n 

A ijk = x j yj 1 j = 1 > 2 »• • • > n 

x k y k 1 k=l,2,...,n, 

iii) Given a finite set of points in the three dimensional space 
p l = ( x l^l,zi,l), ..., P n =( x n»yh> z n>i) j 


yj i 


c k yk 
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the complete invariant system associated with these points is the total 
ity of the four-rowed determinants 


ijkt 


X j 


V; 

z • 

1 


1=1,2,. 

. . ,n 

' 1 
y j 

l 

z j 

1 


j=l,2,- 

. . ,n 

y k 

z k 

1 

y 

k=l,2 t . 

. . , n 


z £ 

1 


£=1,2,. 

. . ,n. 


(!- 6 ) 


All these fundamental invariants are relative invariants of weight one 
under affine transformations. Consider foT example the determinant (F3) 
and the affine transformation 


x' 

y’ 

2 ' 

1 

Then one has 


>1234 = 


a ll a 12 a l 3 

a 14 


X 


X 


a 21 a 22 a 23 

a 24 


y 


y 






= A 


• 

a 31 a 32 a 33 

a 34 


z 


z 


a 4 1 a 42 a 43 

a 44 


1 


1 


x l y{ z i 1 


X 1 y l 

z i 

1 

x* y 2 4 i 

5 

X 2 - 

'2 

Z 2 

1 

x 3 -y$ z 3 1 


X 

y 3 

z 3 

1 

x 4 y\ z 4 1 


X 


^4 

Z 4 

1 


(F7) 


or 

a 1234 = l A l A 1234- (F8 ^ 

The last relation shows that A 12 34 is a relative invariant of weight one. 

In the problem of determinantal loci one has a finite set of points 
in the three dimensional space, namely n station-position points and m 
satellite-position points. Therefore any invariant under affine trans- 
formations, which is a rational integral function of the coordinates of 
these points, must be expressible as a polynomial in the fundamental 
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invariants (F6) . As an example of this take the jacobian matrix for the 
range observations from four ground stations. (see Table E3, Appendix E) 
It has been proved in Appendix C that the determinants of its square sub- 
matrices are invariant under affine transformations. Take a square sub- 
matrix of that matrix, say that of the first twelve rows. That is a 
rational integral function of the coordinates and invariant under affine 
transformations. Therefore it must be expressible as a polynomial in the 
fundamental invariants (F6). This polynomial can be obtained in 
three embedded steps of Laplace expansion: 

1st step - expand the 12-order determinant in Laplace expansion according 
to the first three columns; 

2nd step - expand the 9-order determinant of each term in the 1st step in 
Laplace expansion according to its first three columns, and 
3rd step - expand the 6-order determinant of each term in the 2nd step in 
Laplace expansion according to its first three columns. 

Then one obtains 


A Jl* Jl2 " 2 (-U ai *° 2 * as A] 2 4 1 j 6 A |7 4 J 8 A ^ 9 4 J 1 0 A j 1 1 12 1 . J 2 , ■» 3p j 4 , J 5 , 3 Co J. 7 . 1 


(F9) 


where a^,a 2»°3 result from the above Laplace expansions and P]_>P 2 *P 3 are 
selections of row indices, of these expansions. Aj are given in (E14) of 
Appendix E, and'^’^*^ i s given by 




Z i 


2 


£ 


1 

1 

1 

1 


(F10) 
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If one wants to get an idea about the length of the summation (F9) , 
one will find that there are 


_ r 12 r 9 r 6 _ 12! 9! 6! 

p 1 p 2 P 3 '3 ' 3 ‘ 3 3 1 9 J 3 ! 6 ! 3 ! 3 ! 


369,600 terms. 


The determinant A j equal to zero as it was 

expected, and could be seen by performing the following rank-equivalent 
column-operations in Table E3, of Appendix E: 

i) Add all the x-columns to one of them, say x x -column 

ii) Add all the y-columns to one of them, say y 1 -column 

iii) Add all the z-columns to one of them, say zj-column 
The resulting matrix has its first three columns zero. Thus 


A jl»‘* ■» jl2 


0. 


(EH) 


Such an identity between invariants of the complete system is called a 
syzygy. Syzygies result from the expansions as above, of the determin- 
ants of the submatrices of dimension greater than six, until an admissible 
set of constraints is introduced which leave the jacobian matrix with six 
columns. By constraining the coordinates , yj , z^ ,y2, ^2 > z 3» t ^ ie j ac °bi- an 
for the case of four stations in range networks, (Table 1) is left with 

the columns of x 2 , x 3 >y 3 > x 4 »y 4 j z q * 

Expanding the determinant of its 6x6 submatrix of the first six rows 
as above, one has 


l| ii J« 


PlXp2> 


01+02 




Yl 


M2 T J 2 
Xjj Y Jj 


*4 Y4 *4 1 

X j 4 *j 4 1 

X. Yi Zi 1 

] 5 J S J S 

X)e, Y J 6 l h 1 


4 ' 4 * 4 3 4 4 4 * 4 ‘ 


(F12) 


Notice the presence of fundamental invariants of one dimension: 
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1 


x 2 



and two dimensions: 

i x 3 ^3 1 



Expansion (F12) is of particular importance in the problem of deter- 
minantal loci for range networks. 
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